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Announcing publication of 


The Proceedings of the Symposium on 
—GAS-COOLED REACTORS-— 


held February 10-11, 1960 


sponsored jointly by The Franklin Institute and the Delaware 
Valley Section of the American Nuclear Society. 


Twenty-one papers are featured by well known 
U. S. research people and several foreign scientists 
who have reviewed the developments in this new type 
of atomic power plant. The symposium aimed to 
clarify the problems which must be solved before gas 
cooling can be used to its fullest potential. 


Sir William Cook, C.B., of the U. K. Atomic Energy 
Authority reviewed the British program in gas-cooled 
reactors; a representative from the Swiss firm, Brown, 
Boveri & Cie., discussed their work on a pebble bed 
reactor which they are building for the West German 
Government. Representatives of the Australian Atomic 
Energy Commission and from West Germany pre- 
sented programs being carried out in their respective 
countries. Discussions are included at the end of 
most papers. 


The sessions were devoted to the following topics: 
1 Status of Major Gas-Cooled Reactor Programs 
11 Fuel Elements and Materials 

Il Components and Loop Studies 


Issued as Monograph No. 7 under the auspices of the JOURNAL OF THE FRANK- 
LIN — the Proceedings are now available. Send your order (at $5.00 per 
copy) to 
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20th & Parkway 
Philadelphia 3, Po. 
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FIRST PHONE CALL VIA 
MAN-MADE SATELLITE! 


“Project Echo” sateilite went into a near-perfect circular orbit 1000 miles high, circling 
the earth once every two hours. Its orbital path covered all parts of the U. S. 


BELL TELEPHONE LABORATORIES BOUNCES VOICE OFF SPHERE 
PLACED IN ORBIT A THOUSAND MILES ABOVE THE EARTH 


Think of watching a royal wedding in 
Europe by live TV, or telephoning to 
Singapore or Calcutta—by way of outer- 
space satellites!) A mere dream a few 
years ago, this idea is now a giant step 
closer to reality. 

Bell Telephone Laboratories recently 
took the step by successfully bouncing 
a phone call between its Holmdel, N. J., 
test site and the Jet Propulsion Labora- 
tory of the National Aeronautics and 
Space Administration (NASA ) in Gold- 
stone, California. The reflector was a 
100-foot sphere of aluminized plastic 
orbiting the earth 1000 miles up. 


Dramatic application of telephone science 


Sponsored by NASA, this dramatic ex- 
periment —known as “Project Echo” — 
relied heavily on telephone science for 
its fulfillment . . . 

® The Delta rocket which carried the 
satellite into space was steered into a 
precise orbit by the Bell Laboratories 
Command Guidance System. This is the 
same system which recently guided the 
remarkable Tiros I weather satellite into 
its near-perfect circular orbit 


* To pick up the signals, a_ special 
antenna was used. Previ- 
ously perfected by Bell Laboratories for 
microwave radio relay, it is virtually im- 
mune to common radio “noise” inter- 
ference. The amplifier—also a Labora- 
tories development—was a_ traveling 
wave “maser” with very low noise sus- 
ceptibility. The signals were still further 
protected from noise by a special FM 
receiving technique invented at Bell 
Laboratories. 

“Project Echo” foreshadows the day 
when numerous man-made satellites 
might be in orbit all around the earth, 
acting as 24-hour-a-day relay stations 
for TV programs and phone calls be- 
tween all nations. 


This experiment shows how Bell Lab- 
oratories, as part of the Bell System, is 
working to advance space communica- 
tion. Just as we pioneered in world- 
wide telephone service by radio and 
cable, so we are pioneering now in using 
outer space to improve communica- 
tions on earth. It’s part of our job, and 
we are a long way toward the goal. 


BELL TELEPHONE LABORATORIES 


WORLD CENTER OF COMMUNICATIONS RESEARCH AND DEVELOPMENT 
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SATELLITE ORBITS ABOUT A PLANET WITH 
ROTATIONAL SYMMETRY 


BY 
C. M. PETTY' AND J. V. BREAKWELL' 


ABSTRACT 


A closed-form solution is obtained for the orbit of an Earth satellite which is in- 
fluenced by the third harmonic in the potential in addition to the two dominant ob- 
lateness terms. The results are carried to the ‘‘second order” in the rates of the 
secular variations and in the periodic fluctuations. The secular equations are solved 
in closed form in two cases: when the inclination is close to 63.4° and when it is not. 
In the former case, the stability and period of the slow oscillation of the perigee are 


discussed. 


1. INTRODUCTION 


The objective of this paper is to derive an approximate closed-form 
solution, of high accuracy and without restriction on the inclination 
angle or eccentricity, to the equations of motion of an Earth satellite. 

it has been shown by Roberson (1)? that the motion of a near satellite 
around an oblate Earth may be rather simply described to first order 
in the oblateness by introducing a slowly rotating frame of reference, 
as follows: 


(i) A “mean” orbital plane having constant inclination, rotates 
slowly with constant angular velocity about the polar axis. 

(ii) The angular frame of reference within this plane also rotates 
slowly with constant angular velocity. 


With respect to this doubly moving frame, the motion of the satellite 
differs from pure Keplerian elliptic motion only by perturbations which 
are periodic, including a small periodic motion out of the mean orbital 
plane. 


' Research Scientist and Staff Scientist, respectively, Lockheed Missiles and Space Divi 
sion, Sunnyvale, Calif. 


* The boldface numbers in parentheses refer to the references appended to this paper. 


(Note—The Franklin Institute is not responsible for the statements and opinions advanced by contributors in 
the Journat.) 
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That such a description is possible is, indeed, suggested by more 
classical perturbation procedures, for example, Krause (2), which show 
that the only “‘secular’’ motions due to first-order oblateness are the 
precession of the nodes and the rotation of the line of apsides (at rates 
with which Roberson’s treatment agrees). 

A description somewhat similar to Roberson’s has been developed 
by King-Hele (3), without the small periodic motion out of the reference 
plane. King-Hele’s formulation of the motion of the node, however, 
is unsatisfactory for non-zero “eccentricity.” 

It is the purpose of this paper to obtain a description of the orbit 
beyond first order. Classical procedures indicate that one must expect, 
in general, slow variations in angular momentum, eccentricity, and orbital 
inclination as well as in the line of nodes and the line of apsides. There- 
fore, we examine certain differential equations, derivable from the 
dynamical situation, satisfied by the orbital perturbations as described 
relative to a more general, slowly moving reference frame. As in 
Roberson’s treatment, the slow motion of the reference frame is almost 
completely determined by the requirement that the orbital perturba- 
tions, the solutions of the differential equations, shall be “periodic,” 
that is, all the secular orbit changes are absorbed into the slowly moving 
reference frame. However, in dealing with the interaction between 
the Earth’s “‘pearness” and its oblateness, it is found convenient to 
choose the inclination of the reference plane in such a way that the 
“orbit plane’’ undergoes a slow oscillation about the reference plane 
as the perigee moves around. 

We will assume that the Earth has rotational symmetry and that 
the motion of a satellite is governed solely by the Earth’s potential 
function; that is, other forces like drag and attraction of the sun and 
moon are not considered. We introduce a right-handed inertial rec- 
tangular coordinate system with the origin at the center of mass of the 
Earth and the z-axis coincident with the rotational axis and positive in 
the north direction. The positive x-axis is chosen in the direction of 
the vernal equinox (of, say, 1950). The potential V(r, 6) of a point 
in space is a function of the distance r to the center of mass and de- 
clination 6. The function V is given by 


4 
(3 30s;? + 35s,*) ee | (1) 


where R, is the equatorial radius of the Earth and yu, Js, Js, J4 are con- 
stants. We abbreviate sin 6 by s;. The odd harmonic involving J; is 
due to the possible lack of symmetry of the Earth with respect to the 
equatorial plane. 
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In deriving a closed-form solution, the accuracy and complexity of 
the results and the amount of effort required depend in a large measure 
on the choice of dependent variables and the choice of an auxiliary inde- 
pendent variable. These variables, in turn, should determine the 
position of the satellite in terms of some commonly accepted and 
readily understood parameters such as right ascension \, declination 6, 
and radial distance r to the center of the Earth. The end result then 
should be judged on the complexity and accuracy with which X, 6, and r 
are determined. 

On the other hand, the constants which appear in a closed form 
solution do not have to be simply related to initial conditions if the 
purpose of the solution is for tracking and prediction of a satellite. In 
this case they are simply constants to be determined by observations 
using some statistical technique. 


Il. EQUATIONS OF MOTION 


The unit vectors along the positive x-, y-, and z-axes are denoted 
by i, j, and k, respectively. A moving right-handed orthogonal set of 
unit vectors is introduced by 


where Q is the Euler matrix 


— SeSaCy + SuCeCy 


— — — + CuCaCy 


— 


with Euler angles Q, y, and w. 

It may be helpful for the reader to visualize 2, y, w as prescribed 
functions of time such that @ is the right ascension of the ascending 
node, y the inclination angle, and w the argument of perigee. However, 
these functions are chosen in a very special way, and their exact de- 
termination is not obtained until the end of this paper. Actually, they 
represent these (vaguely defined) physical parameters only in some 
sort of average sense. A second moving set of right-handed orthogonal 
unit vectors e,*, e.*, e,* is introduced by 


. 
1 
: 
(3) 
: 
| = | e: | (4) 
| @3 { @3 | : 
: 
as 
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where 


— — 


in the transformation matrix to spherical coordinates relative to the 
set 

Again it may be helpful for the reader to visualize @ as the true 
anomaly and ¢ as the angular deviation, as seen from the origin, of the 
satellite from the slowly moving plane determined by e, and e». 

We will choose @ as the auxiliary independent variable and, by 

primes, denote derivatives with respect to 6, for example, w’ = 48 A 
dot over a dependent variable will denote its derivative with respect to 
dr 
dt 
The variables @ and ¢@ are determined such that 


time, for example, * = 


r= re,* 
where r is the radius vector of the satellite and 
Ae,* + Be.* + Ce;*. 


Since e,* are orthogonal unit vectors, 


where 


| — | 


| 
| 
| 


Krom Eqs. 7 and 8 we have 


2ia* + ra* — B*5*r = B 


+ rp* + = C. 
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Q* = | — Se Ce 0) | (5) 
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e; a B | 
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a* p* | 
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The three dependent variables which will be chosen are p, ¢, and + 
given by 


(11) 


One can then compute 


tp*(1 + + F; + Fio) | 
= + F,] | 
rp*(1 + + F;)o — 
where 
Fi = — 


= ¥' — 


| 


F; w’ + 


The differential equations (10) may now be written in terms of the de- 
pendent variables p, ¢ and r. 


pT 
T 


p + p(l + + 


+ p + 2F.0' 4s Fi Fw + al F, de F.’) 4 F yr 


T 


(aC — A)(1 + 
(16) 


a” +oa(1 + F;)? + F, F + (a’ + 


+ 2F.,00' — + + F;) 


(17) 
T 
(1 + F3)r’ + Filor)’ + r(F;' + F,F:) 


B 
+ + ra( fF,’ — F, — = (18) 
p*r 


The above are the general equations of motion in our coordinate 


p= 1 
o = tang | 
T= bp 
( 
12 ) 
(13) 
) : 
| 
| 
5) 
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system but we will consider only the case corresponding to the Earth’s 
potential, Eq. 1, where 


(r, 4) 


’ 


or 


V(r, 6) 
06 


TC; 


(1 + g?) — 6) 
06 


TCs 


The radial distance 7, right ascension \, and declination 6 of the 
satellite are given by 


r = + 


Ch Le 
ss = (1 + + 


Up to this point no approximations have been made and the differen- 
tial equations and other expressions are exact. In the next section 


Eqs. 16, 17, and 18 will be solved approximately when A, B, and C 
are given by Eq. 19. 


Ill. A CLOSED-FORM SOLUTION 


If the potential V is given simply by V = uw/r and Q, y, and w are 
chosen to be constant, then Eqs. 16, 17, and 18 are given by 


2)—3/2 
(1 + o?)—3/2, 


We may then choose as a closed-form solution to Eqs. 16, 17, 
and 18 


p = — (1 + &»), 


(22) 


{ | 
i = ' 
| 
= 19) a 
“+ 
| 
is 
} 
= 
” (21 ) 
+o = 0, | 
og =). 
} 
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where & = h’ = 0 and by Eq. 20 


=p! 


2+ tan | | 
Case j 
| 


The particular constants Q, y, w, €, # which enable any solution to be 
written in this form are then to be determined by observation of the 
satellite. 

Returning to the equations with V given by Eq. 1, the problem is 
to determine 2, y, w, €, # as slowly varying functions of @ such that any 
solution to Eqs. 16, 17, and 18 has the form 
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where aj, dj, 6, 6;, d;, d; in the expressions for p;, ¢; and 7; are expressed 
in ascending powers of the J,’s with coefficients which are functions of 
2, y, w, e«, and h. The four terms ay, d;, 5; and 6; play a very special 
role which we will discuss presently. 

We adopt the following convention as to which terms of the solu- 
tion are to be retained. Terms with J; as a factor will be considered 
as first order, those with J,*, J;, or J, as a factor will be considered 
second order. Terms higher than the second order will not be retained 
in this paper, consequently V in Eq. 1 will be restricted to those terms 
explicitly written. For conciseness, we will write 


(27) 


It is clear that Q’, y’, w’, ¢, h’ and aj, a;, b;, 6;, d;, d; are of the first 
order, at least. For the moment, let us consider terms of, at most, the 
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first order. The equations then become 


Replacing p, ¢, tr by Eqs. 24, 25 and 26 and equating coefficients, we 
find that 2’, w’, di, ds, ds, bo, Go, G2, G2, As, Gs, Gy are the only terms 
of the first order and that among the derivatives of these terms, only 0,’, 
Ge’, a3’, a4’, are of the second order, the others being at 
least the third order. 

Proceeding now to a second-order analysis, the equations become 
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It is now necessary to discuss the basic philosophy which motivates 
our approach to this problem. For the moment, suppose a, di, b;, 6, 
are set equal to zero in Eqs. 24 and 25. Our procedure then would be 
to choose the rates 2’, y’, w’, €’, #’ such that they eliminate the resonance 
(first harmonic terms in @) in Eqs. 16 and 17 and the constant (slowly 
varying) term in Eq. 18. In this manner we obtain five ‘‘secular”’ 
equations which may strongly interact with each other for critical 
values of the variables (for example, y near 63.4°), but the functions 
pi, 71, 7, Will in all cases have the same form. The problem is then 
completely transferred to the solution of these secular equations. 
However, a basic difficulty arises from this procedure which, as far 
as we have carried it, is due entirely to the presence of the odd harmonic 
term involving J; in the potential function. 

When we gather the coefficient of se, in Eq. 17, the dominant term 
in Q’ is 2s,Q’, but also there will appear terms involving J; which do 
not have s, as a factor. Similarly, in Eq. 16 the coefficient of cy will 
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have the dominant term 2e(w’ + c,Q’) and_also certain J; terms which 
do not have ¢ as a factor. This procedure then leads to secular equa- 
tions which strongly interact for small eccentricities or small inclination 
angles. Moreover, the unfortunate feature is that in the calculation 
of the coefficients in p;, ¢:, and 7; as explicit functions of ¢€, w, y, and h 
there will appear terms which have either ¢ or s, in the denominator. 
To the second order, this does not occur, but these would appear if 
third-order terms had been retained. However, these higher order 
terms should also be retained when either € or y is small. Consequently, 
the treatment is unsatisfactory in these two cases. 

This difficulty is removed by introducing the terms d,s, in Eq. 24 
and 6\ce,. in Eq. 25. The derivative 6,’ will appear in the coefficient 
of ses. in Eq. 17. Its purpose is to remove all terms which do not have 
s, asa factor. Similarly, the derivative 4,’ will appear in the coefficient 
of cs in Eq. 16, and its purpose is to remove all terms which do not have 
€ as a factor. 

The additional introduction of aycs in Eq. 24 and bse... in Eq. 25 is 
not an essential feature, but it is convenient to do so. For, if we did 
not do this, we would find that the amplitude of the fluctuation in € or 
does not become small when the corresponding variable is small, again 
due to the presence of the J; terms. This means, for instance, that for 
small eccentricity it would, at times, be necessary to interpret € as 
being negative. Therefore, a; and }; are used to keep the rate e’ pro- 
portional to e and y’ proportional to s,. 

Summarizing the philosophy of our procedure, we require that ay, 
4, b;, 6; be functions of y, w, e, and h such that every remaining coeffi- 
cient of a first harmonic term in @ in Eq. 16 has ¢ as a factor, and every 
remaining coefficient of a first harmonic term in Eq. 17 has s, as a factor. 
The rates w’, 2’, e’, y’, and h’ are then chosen such that they remove 
all remaining first harmonic terms in @ in Eqs. 16 and 17 and also the 
constant (slowly varying) terms in Eq. 18. 

Now the choice of ai, d:, 6;, 6: is not necessarily unique, for it is 
possible for instance, to absorb into these variables some additional 
terms in Eq. 17 having s, as a factor. The choice that was made is 
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and consequently 


, 3; 
+ = — ~>8,(1 — 
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Returning now to the second-order equations, we first solve Eq. 31, 
then Eq. 30 and finally Eq. 29. In the computation, },, bi, ai, a; are 
treated as first-order terms. Consequently, the reader will find them 
involved in the computation of d,, di, ds, d2, ds, ds, bo, b2, Go, 
3, G3, a4, Gy in the J; terms. In the following pages of our results we 
have evaluated the above coefficients explicitly according to our choice 
in Eq. 32. 

The “secular equations” are given by: 
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The coefficients d,, d, in Eq. 26 are given by: 
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The coefficients 4;, 6, in Eq. 25 are given by: 
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In the next section attention will be directed toward the secular 
equations (34 to 38). 


IV. SOLUTION OF THE SECULAR EQUATIONS 


These non-linear equations demonstrate an interesting interaction 
when the inclination angle y is in the neighborhood of 63.4°. In this 
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case, the term /.(2 — 


2 
havior of w, and there is a strong interaction with Eq. 36. If we write 
y = tan-'2 + J. and ignore terms of the third order or greater in 
Eqs. 35, 37 and 38 we obtain 


s*,) in Eq. 35 no longer dominates the be- 
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and if J's not close to 18’ the dominant term for n’ is given by 


30 3572) (82) 


The only possible equilibrium positions of this slow motion for 
0 corresponds to s». = 0. In this case e, and are constants, 
where n has the precise value which makes w’ = 0. Now Eqs. 80 and 
81 may be combined with Eq. 82 and ¢ and h solved exactly as functions 
of 7. These solutions justify, to our order of approximation, treating 
e and hk as constants in the solution of Eq. 79 combined with Eq. 8&2. 
Let 
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Phe two cases of equilibrium correspond to w = + 5 and w = 0, x. 


The stability of the equilibrium positions is decided by linearizing 
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Eq. 83 in the neighborhood of an equilibrium point. We find that 
= + ~are stable equilibrium positions if J,/J2 > 35/18, and w = 0, 


2 
x are stable equilibrium positions if J,/J. < 35/18. The period of the 
small oscillations is very long, being 


V30 
35 


orbital revolutions. 
Returning to Eq. 83, we may solve for w as a function of 7 by solving 
the linear equation 
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where A is the constant of integration. A family of curves of w versus 


Ay ( R ) is given in Fig. 1. In this case, we have chosen J, = 1.637 
uk, 


35 h.? 
x 10-4, J, = 1 10-5 (> + se) and (45) = 0.05, and neg- 


lected the terms in a, band c. With this neglect it should be noted 
that for the family of curves, Eq. 86, with parameter K involves the 
additional parameter ¢,/,? but not ¢, and h, separately. 

The above phenomenon is further discussed in (4) where the inter- 
action with “sideways drag”’ is also considered. It should be noted 
that the Azin (4) wasi — tan~' 2, with signifying the average (wrt @) of 
the instantaneous inclination angles. This average 7 is related to our 
present parameter y by: 


The term involving J;/./. shifts the northern equilibrium position rela- 
tive to the southern one (4, p. 4). Figure 1 of (4), however, was com- 
puted with J; = 0 and thus differs from our present Fig. 1 only by a 
shift in the origin of abscissae. 

When ¥ is not close to 63.4°, w is a monotone function, and we may 
rewrite Eqs. 36, 37 and 38 with w as the independent variable. These 
three equations may now be solved approximately by treating y, «and h 
as constants on the right-hand side. We then obtain: 
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where y., h., € are constants and J; = J;(R.y/h.2)'. We are not 


justified in carrying terms higher than the first order in Eqs. 87, 88 and 
89 since Eqs. 34 through 38 were carried only to the second order. 

We may now substitute Eqs. 87 through 89 in Eqs. 34 and 35 and, 
dropping terms higher than the second order, these may be integrated. 
We then obtain 
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where w, and Q, are constants. 

The five constants y,, €, /., @., 2. carry the burden of enabling the 
total solution (these secular equations adjoined to the results of the 
previous sections) to be expressed in the form which we have adopted. 
It should be observed that, in general, the functions y, ¢, 4, w and Q 
do not have the values y., €o, Ao, @., 2, when 6 = 0, although the form 
of Eqs. 79 through 8&3 could have been chosen so that this is the 
case. Similarly, in the previous section the functions a;, di, }1, 6: 
were solutions of certain “‘secular’’ differential equations, and their 
form could have been chosen such that they assume the value zero 
when 6 = 0. We have not done this simply because it would have 
complicated the appearance of a number of other terms (for example, 
d,, dy, ds, etc.). As far as the authors can see, this procedure gains 
nothing. 

In order to obtain a complete second-order theory, it is clear that it 
is necessary to carry the secular differential equations to the third order. 

It is interesting to interpret geometrically the major effects of the 
J;terms. It has been estimated by Kozai (5) that 


for the planet Earth. The negative value of J; implies a pear-shaped 
Earth with a bulge in the southern hemisphere. The plane P de- 
termined by Q2 and y¥ will always contain, as it rotates, the major portion 
of this bulge on the same side, provided the orbit is not polar. Now 
the presence of the term }(/;/J2)c,ecs in the expression for ¢ implies 
that the orbit is tilted with respect to P so that the perigee remains on 
the side toward the bulge and the apogee remains on the opposite side. 
This tilt then follows the perigee as it precesses. An observer who 
calculates a best-fitting plane to the orbit would interpret this as a 
sinusoidal fluctuation of the inclination angle of period 27 in w. 

The term —3(J;/J2)s,(1 — &)se,. in the expression for p, shifts 
the orbit in the direction of the north-south apexes, so that the satellite 
comes closer to the center of mass in the northern hemisphere and farther 


J 3 
2 
- 
/ J, (6 
+ - + 5 + 3 s?,,) |} sw (91) 
3 
: 
ig 
ue 
6 
J; 2.2 x 10 (92) 
: 
; 
; 


Oct., 1960.] SATELLITE ORBITS 281 


away in the southern hemisphere (the one containing the bulge). If 
an observer calculates a best-fitting ellipse to the orbit, then as the 
line of apsides rotates he would interpret this effect as a sinusoidal 
fluctuation of the eccentricity of period 27 in w such that the perigee 
distance is least when it is at the northern apex. 

There is a small but interesting effect due to a north-south balance 
of forces. Intuitively, due to the rotational symmetry, it is clear that 
we can have a circular orbit with center on the rotational axis and lying 
in a plane perpendicular to this axis, provided the latitude of the plane 
is adjusted so that the north-south forces just balance for points on the 
orbit. If we set ¢,=y.=0 in our equations, then we obtain ¢ = (3/2) J;. 
Using the value of J; given by Eq. 92, it is found that a near-Earth 
orbit of this type would remain about 70 feet south of the equatorial 
plane. In this paper the equatorial plane has been defined as that 
plane through the center of mass perpendicular to the rotational axis. 
An observer on the ground whose plumb bob is perpendicular to the 
rotational axis would find this orbit slightly to his north. 

If J; is second order compared with J, as we have assumed in this 
paper, then we should be able to approximate the solution up to the 
second order for at least several revolutions without the J; terms. 
For example, suppose we set J; = 0 everywhere that it occurs in our 
solution and we let 2,*, €.*, 4.*, y.* and w,* be the five constants in 


Eqs. 87-91. Also denote the variable @ by 6* in this solution. We will 
suppose for the moment that ¢,, s,, are not too small, and we choose 
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solutions agree up to the second order in radial distance, right ascension, 
and declination. 

When either the eccentricity or s, is small, the representation of the 
motion without the J; terms can be troublesome. Suppose, for ex- 
ample, we consider a polar orbit with e, = 0. In this case, our solution 
gives an orbit which closes after one revolution. Physically, there is a 
balance in the gravitational attraction such that an orbit which cuts 
the north polar axis at right angles also cuts the south polar axis at 
right angles. Without the J; terms, a good approximation is obtained 
by choosing the perigee on the north polar axis and ¢,* = (—3/4)(J3/J:). 
However, the rotation of the line of apsides causes the perigee to move 
away, but the situation is actually static. 

A similar case arises when y, = 0 but €,is rather large. The perigee, 
in this case, remains at the southern apex of the orbit as the perigee 
precesses in right ascension. Without the J; terms, a good approxima- 
tion is obtained by choosing the inclination angle y,* = (—3/4)(J3/Js)e, 
and the perigee at the southern apex. However, the rotation of the 
line of apsides causes the perigee to move toward the northern hemi- 
sphere when actually the perigee remains at a constant southern latitude. 
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ELECTRIFICATION OF CARBON BLACK BY CONTACT 
WITH A METAL SURFACE 


BY 
CHARLES J. BRASEFIELD' 


ABSTRACT 


A measured quantity of carbon black was allowed to slide down a steep surface of 
nickel (presumably oxidized) which was electrically grounded and the charge ac- 
quired by the carbon black was determined. Blacks were selected whose oxygen 
content varied from 0.4 to 11.6 per cent. It was found that a black of low oxygen 
content acquired a large positive charge and that the charge decreased with increasing 
oxygen content, becoming zero for 4 per cent oxygen and negative for larger oxygen 
contents. It thus appears that the electron affinity of carbon black increases with 
increasing oxygen content of the black. The positive charge acquired by a black of 
low oxygen content was found to decrease as the black was heated to higher tempera- 
tures, presumably because of increasing oxygen content of the black. The oxidation 
of certain blacks is not permanent if the temperature reaches only 400° C. but is 
permanent if it reaches 600°C. The oxygen content at a particular temperature 
ceases to increase after the first 30 seconds at that temperature. The oxidized surface 
of carbon black which is responsible for its high electron affinity cannot be produced 
by prolonged exposure of a black to atmospheres of CO, COs, or CO>. 

Cigarette smoke is usually negatively charged, but if the cigarette is heated to 
300° C., the smoke is positively charged, as though the electron affinity of the tobacco 
had been increased by oxidation. 


1. INTRODUCTION 


It has been shown (1)? that carbon particles in the exhaust of a motor 
vehicle are negatively charged. This charge was detected by means of a 
cylindrical Faraday cage of window screening, 50 cm. in diameter and 
50 cm.long. In the center of the cage was a polonium button connected 
directly to the input terminal of a Keithley electrometer (Model 200A, 
zero centered) on which the cage rested. When positive (or negative) 
space charge enters the cage, the potential of the button may be several 
volts positive (or negative) with respect to ground. During the course 
of the experiments with motor vehicles, cigarette smoke was frequently 
blown into the Faraday cage. The electrometer deflection was some- 
times negative and sometimes positive, but if short fast ‘“‘drags’’ were 
taken, the cigarette smoke was generally negatively charged. 

In order to learn more about the process by which space charge is 
generated by motor vehicles and by cigarettes, it was decided to in- 
vestigate the charge acquired by carbon black after sliding over a 
grounded metal surface. Thomas (2) has shown that a stream of fine 
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coal particles acquires a positive charge after contact with a metal 
surface, but if the coal is oxidized by heating in air at 350° C. for some 
minutes, the charge acquired is negative and its magnitude depends on 
the degree of oxidation of the coal. One purpose of the experiments 
planned with carbon black was to investigate the relationship between 
the oxygen content of the carbon black and the charge acquired, and, 
if possible, to learn something about the nature of the oxidized surface 
layer. 


THERMOCOUPLE 


HEATING COILS 


—+-METAL DISCS 


+-BAFFLE 


+—-NICKEL 


O1294:5 


_|*—ELECTROMETER CUP 


Fic. 1. Furnace, with ten metal discs inserted 


2. APPARATUS AND PROCEDURE 


The experiments were carried out under conditions of pressure and 
humidity that existed in the laboratory. The procedure consisted 
simply in allowing a measured quantity of carbon black to slide down a 
steep surface of nickel (presumably oxidized) which was electrically 
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grounded. ‘The carbon black was then collected by a shallow metal 
cup which rested on the input terminal of a Keithley electrometer 
(model 210, scale 4-0-4 volts). The electrometer deflection was pro- 
portional to the charge acquired by the carbon black. Generally, it 
was desired to heat-treat the carbon black before it was allowed to 
slide down the nickel surface. Several ways of doing this were tried 
and the arrangement illustrated in Fig. 1 was finally adopted. The 
furnace and electrometer cup were enclosed in a cage of wire screening 
which was grounded. 

The lower end of the furnace was conical on the inside and was 
lined with nickel. A nickel baffle was inserted to prevent any carbon 
black falling through the cone without striking the walls. Samples of 
carbon black were placed on each of seven metal discs. The tempera- 
ture of the upper three discs was considerably below average, so they 
were not used. The temperatures of the remaining seven discs were 
nearly the same, being slightly higher than the temperature indicated 
by the thermocouple. The discs were supported by a hinge on one side 
and by a pin in a vertical rod on the other side, so that ten pins were 
required to support the ten discs. These pins were inserted in the rod 
at different azimuths. If the rod was turned about 15°, the lowest 
disc was free to fall. Another turn of 15° would release the next disc, 
and soon. In this way, seven samples of carbon black could be made 
to slide down the nickel cone. The volume of carbon black in each 
sample was the same, called a standard dose. It was the amount held 
by a spoon 16 mm. in diameter and 3 mm. deep. The bottom of the 
spoon was made of window screening so that the carbon black could be 
dusted uniformly over a disc. 

On a given day, the charges acquired by successive standard doses 
of the same black generally differed by less than 20 per cent. When- 
ever possible, the average of at least three readings was taken. The 
charge tended to be less when the relative humidity was high. The 
data here reported were taken on days of low to moderate humidity 
(40-60 per cent). 


3. RESULTS 
(a) Vartation of Charge with Oxygen Content of Carbon Black 


In order to investigate the relationship between the oxygen content 
of the carbon black and the charge acquired, a number of blacks were 
selected? whose oxygen content was known (3). <A different black was 
placed on each of the seven discs, except one black was repeated. The 
charge acquired when a standard dose of the black was allowed to slide 
over the (oxidized) nickel surface at room temperature was measured 
three times on the same day. The mean charge is shown in Fig. 2 as a 


* Kindly furnished by George E. Marsh, Manager, Special Blacks Division, Godfrey L. 
Cabot, Inc. 


é 
: 
: 
| 
‘ 
08 
: = 


286 CHARLES J. BRASEFIELD 


function of the oxygen content of the various blacks. ‘There is a five- 
fold variation in the apparent density of these blacks, Carbolac 1 being 
the least dense. If account is taken of this variation in density, the 
charge acquired by unit mass of carbon black can be computed. This 
is indicated by the dashed curve of Fig. 2. It appears that the charge 
becomes more negative as the oxygen content of the carbon black 
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hic. 2, Variation, with oxygen content of the black, of mean charge acquired by a 
standard dose of carbon black on sliding over an (oxidized) nickel surface The charge ac- 


quired per unit mass of carbon black is indicated by the dashed curve 


increases. This is what would be expected if the work function (elec- 
tron affinity) of the carbon increased with increasing oxygen content, 
if the same way that oxidation of nickel increases its work function (4). 

If the surface over which the carbon black slides is a metal other 
than nickel, with a work function lower than that of nickel, the same 
sort of curve is obtained except that the horizontal line indicating zero 
electrometer reading is moved upward. For example, after sliding down 
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an aluminum chute, the charges acquired by unit mass of Gastex, 
Elf 0 and Carbolac 1 are —2.0, —2.4 and —4.2, respectively. This 
indicates that the horizontal line corresponding to zero electrometer 
reading for aluminum is to be found on Fig. 2 at an ordinate of 
about +2.6. 


(b) Variation of Charge with Particle Size 


Since the particle size generally decreases with increasing oxygen 
content, it might be argued that the independent variable in Fig. 2 
should be particle size rather than oxygen content. This is not so, for 
several reasons. The particle sizes of Elf 4, Elf 0 and Mogul A are 
practically identical, but there is a marked difference in the charge 
acquired by a standard dose. Furthermore, Supercarbovar and 
Carbolac 46 have about equal particle diameters (about half that of 
Elf 0) but the mean charge for Supercarbovar (low in oxygen) was +0.7 
while the charge for Carbolac 46 (high oxygen content) was —0.6. 
Elftex 5 has about the same oxygen content as Gastex but the particle 
diameter of Elftex is less than half that of Gastex. Nevertheless, the 
charges acquired by unit mass of Elftex and Gastex are essentially the 
same. It must be concluded that if particle size has any effect on the 
charge acquired by a standard dose of carbon black, the effect is of 
secondary importance. 


(c) Variation of Charge with Heat Treatment 


In order to investigate the nature of the oxidized surface layer, 
experiments were undertaken to see whether the oxygen content of a 
carbon black could be increased by raising the temperature of the black 
for certain periods of time in the same way that Thomas (2) was able to 
increase the oxygen content of fine coal. Standard doses of Gastex 
were placed on the seven discs, the furnace heated to a maximum tem- 
perature of 400° C., and then allowed to cool. The charge acquired 
when the Gastex was dumped at various temperatures, both when the 
temperature was increasing and when it was decreasing, is shown in 
Fig. 3. It appears that the electron affinity (and presumably the 
oxygen content) of Gastex increases with increasing temperature, 
becoming nearly constant at 400° C., at which temperature the electron 
affinity of the (oxidized) Gastex is identical with that of the (oxidized) 
nickel surface over which the Gastex. slid. However, the oxygen 
content of the Gastex is not permanently increased, since the charge 
acquired by a standard dose at 70° and at 50° is about the same for 
falling temperatures as for rising temperatures. If, on the other hand, 


the heating is carried to 609° C., then the oxidation of the Gastex ap- 


* Bosworth (4) reports that the work function of nickel is 5.0 ev. and that of nickel oxide 
is 6.4 ey 
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pears to be permanent (4-5 per cent); in this case, on cooling to room 
temperature, the charge acquired is slightly negative. Sterling MT is 
also permanently oxidized after heating to 600° C, 

Figure 3 is interesting for another reason. It might be argued that 
the change in electron affinity of the carbon black with rising tempera- 
ture was due to the loss of a surface layer of water molecules, rather 
than to an increase in oxygen content. If this were the case, however, 
the curve for falling temperatures should be quite different from the 
curve for rising temperatures. Furthermore, tests show that thin 
layers of carbon black can be stored for 70 hr. in a desiccator over con- 
centrated H.SO, without appreciably changing the charge acquired 
by a standard dose. 
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Fic. 3. Variation with temperature of the charge acquired by a standard dose of Gastex, both 
for increasing (0) and for decreasing (x) temperatures. 


Thomas (2) found that the oxidation of coal at 350° C. increased 
continuously over a 15-min. period. An experiment was therefore 
carried out to see whether the oxidation of carbon black would increase 
with prolonged heating. The furnace shown in Fig. 1 was heated to 
400°C. The top of the furnace was opened long enough to insert 
seven discs loaded with Gastex. The lowest disc was dropped about 
30 seconds after insertion in the furnace, and the electrometer read 
—(.1 division. The remaining discs were dropped at 15-min. intervals, 
and theelectrometer readings were +0.2,0, —0.3, —0.1, —O0.15 and —0.1, 
respectively. Thus the reading after 90 minutes is the same as the 
reading after 30 seconds, indicating that in this case, the oxidation of the 
carbon black was complete within 30 seconds. These samples of carbon 
black were saved and tested again the next day. The electrometer 
readings were the same as for fresh Gastex, indicating that the oxidation 
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at 400° C. was not permanent. However, a sample of Gastex which 
had been heated to 600° C. produced a negative deflection of the elec- 
trometer when tested the next day, showing permanent oxidation. 

Measurements were made of the charge acquired by samples of Elf 0 
at various temperatures up to 600°C. In all cases the charge was 
essentially the same, suggesting that the oxygen content of Elf 0 
(5 per cent) is not increased by heating to 600° C. 


(d) Experiments with CO, CO, and C,O, 


An attempt was made to produce an oxidized surface layer on carbon 
black by soaking thin layers of the black for 24 hr. in an atmosphere of 
(a)® CO, (b) COs, and (c)** C,O,. However, the charge acquired by 
carbon black thus treated was not appreciably different from the 
charge of untreated carbon black. 


4. GENERATION OF SPACE CHARGE BY MOTOR VEHICLES 


To simulate conditions in the exhaust pipe of an automobile which 
is operating with its engine choked, Sterling MT was made to slide over 
a piece of rusty iron. Sterling MT is a thermal black which, because 
of its method of production, should be very similar to the carbon from 
a choked engine. A standard dose of Sterling MT, after contact with 
rusty iron, acquired an average charge of +0.15 unit. It therefore 
appears that the carbon from a choked engine is somewhat more highly 
oxidized than Sterling MT, since tests show the former has a negative 
charge (1). When the automobile engine is operated under normal 
conditions, any carbon in the exhaust should be still more oxidized and 
therefore have a greater negative charge. Since under normal condi- 
tions the exhaust is generally positive (1), it is reasonable to ascribe this 
charge to water droplets. 


5. GENERATION OF SPACE CHARGE BY CIGARETTES 


It was noted above that cigarette smoke is generally negatively 
charged. This can be explained by assuming that the electron affinity 
of the carbon particles in the smoke is greater than that of the tobacco 
over which the smoke passes. If the cigarette is heated to 300° C., 
either before or while it is being consumed, the smoke is generally 
positively charged, as if the electron affinity of the tobacco had been 
increased by oxidation. The loss of water by the cigarette on heating 
cannot explain the effect, because a cigarette which has been in a 
desiccator over concentrated H.SO, for 20 hr. still produces negatively 
charged smoke. 


5 Kindly prepared by R. F. Trimble of the Chemistry Department of Southern Illinois 


University 
* Kindly prepared by K. A. Van Lente of the Chemistry Department of Southern Illinois 


University. 
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6. CONCLUSIONS 


The charge acquired by a measured quantity of carbon black, after 
sliding over a nickel surface, becomes more negative as the oxygen 
content of the black increases, presumably because the electron affinity 
of the carbon black increases with increasing oxygen content. The 
effects observed cannot be accounted for either by differences in particle 
size of the blacks or by surface layers of water molecules. The oxygen 
content of a black can be increased up to a total of about 5 per cent by 
heating the black to 400—600° C. The oxidation may not be permanent 
if the temperature reaches only 400° C. Prolonged heating at a par- 
ticular temperature does not produce any more oxidation than is pro- 
duced in the first 30 seconds at that temperature. The oxidized surface 
of carbon black which is responsible for its high electron affinity cannot 
be produced by prolonged exposure of a black to atmospheres of CQO, 
CO, or C;,O.. If a cigarette is heated to 300° C., the smoke is positive 
instead of negative, as if the electron affinity of the tobacco had been 
increased by oxidation. 
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ON THE ASSOCIATED FLOW RULE OF PLASTICITY 
BASED ON CRYSTAL SLIPS* 


BY 


T.. 


ABSTRACT 


Associated flow rule has been commonly assumed in the limit analysis of metal 


structures. In this paper, aggregates of differently orientated metal crystals with finite 
elastic constants are considered. The coincidence of the plastic potential with yield 
surface of these polycrystalline aggregates is theoretically derived from the slip 
characteristics of single crystals. The deviation of the plastic potential from the 


yield surface of the aggregate increases with the elastic anisotropy of the individual 


crystals 


INTRODUCTION 


To simplify the calculation of the ultimate load that the structure 
can carry, simple yield surface and its associated flow rule are commonly 
assumed. The present paper does not attempt to derive the yield sur- 
face under different histories of loading but aims to show the condition 
required of the individual crystals in the aggregate to give the coin- 
cidence of plastic potential with yield surface. 

The state of stress at a material point may be represented by six 
independent stress components 7,; of a stress tensor. Using these six 
stress components as coordinates in a six dimensional space, a state of 
stress can be represented by a point in this stress space. The region 
in the stress space, within which only elastic changes in strain occur, 
is bounded by a surface called yield surface. For work-hardening 
materials, when the stress point goes beyond the existing yield surface, 
both elastic and plastic strains are produced and a new yield surface 
F is formed. <A similar strain space can be constructed by using the six 


independent strain components ; as coordinates, e;; = 
OX ; Ox, 


where x,, x, and x; are a set of rectangular coordinates and 1, is dis- 
placement along x,-axis. The strain coordinates are superposed on the 
corresponding stress coordinates. In order for the dot product of the 
stress vector and the incremental strain vector to give incremental work, 
C11, Cor, Cas, 2€12, 22, and 2e,; are used as strain coordinates. A hypoth- 


* This research was supported by the United States Air Force through the Solid State 
Science Branch of the Office of Scientific Research of the Air Research and Development 
Command under Contract AF 49(638)-20. 

Professor of Engineering, University of California, Los Angeles, Calif. 
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esis in plasticity states that there exists a function of stresses g(7,,) such 
that the plastic strain increments are given by 


Og 
OT, j 


;, dy 


where dd is a non-negative infinitesimal and double prime denotes 
plastic. This equation is called flow rule. A further hypothesis is 
made that the yield surface f and the potential function g be identical. 
This is called the associated flow rule of plasticity. It has been shown 
by Hill (1),? Koiter (2) and Drucker (3,4) that these hypotheses give 
uniqueness of solution of stress and strain field for elastic plastic bodies. 
Drucker (3) has shown that this associated flow rule is satisfied for 
materials with positive plastic work done by the incremental stress 
dr,de"’,,; > 0. This property dr,,de’;; > 0 does not follow directly 
from the present conventional laws of mechanics or thermodynamics as 
pointed out by Hill (5), Bland (6) and Drucker (4) and hence is not 
common to all materials. Since metallic structures are commonly 
used, it is the aim of the present paper to study the validity of this 
associated flow rule for polycrystalline metals. 

Single crystal tests by Taylor and Elam (7, 8) have shown that 
under high stress, slips occur along crystallographic directions on 
crystal planes. These slips are the main mechanism of plastic strain. 
For face-centered cubic crystals such as aluminum and copper there are 
four crystal planes and on each plane there are three slip directions at 
60 deg. apart, so there are twelve slip systems. It has been observed 
that single slip depends on the resolved shear stress (shear stress along 
the direction of slip on the crystal plane) and is independent of the 
normal pressure on the plane. ‘The resolved shear stress to cause fur- 
ther slip is termed the critical shear stress. The slip in one system 
generally not only gives strain hardening, that is, raising the critical 
shear stress in the particular active slip system, but also gives hardening 
to the inactive slip systems. This is known as latent hardening. In 
most metals there are both active and latent hardenings. The relative 
amounts of these hardenings vary for different metals. 


Associated Flow Rule of Plasticity for Single Crystals 


Let 7., be the critical shear stress of the slip system » in the crystal, 
de’’,;, a given plastic strain increment and @,;, a stress which will produce 
this strain and o*,; be any other stress within the yield surface. Let 
r, and r*, be the resolved shear stresses of the active slip systems corre- 
sponding to o,; and o*,,;, respectively. 7,’s in the active slip systems 
are equal to the critical. Since o*,; is within the yield surface, r*, is 


? The boldface numbers in parentheses refer to the references appended to this paper. 
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always less than the critical. Hence 


(r, — 7*,)dy". > 0. 


The incremental work due to plastic strain increment 


dw — dw* = — o*,;)de”’;; = ¥ (1, — > 0. 
nel 


to 
wate 


Writing (¢,, — o*,,) as do,;, the above equation gives do,de’’,; > 0. 
This is similar to Drucker’s condition for stable material (9). (¢,;—@,,*) 
makes an acute angle with de”’,;. Hence de’’,; is normal to the yield 
surface for single crystals. This has been given by Bishop and Hill (10) 
as the principle of maximum work which shows the coincidence of 
plastic potential with yield surface for single crystals. 

Let y:, v2 and y,; be the rectangular coordinates with one axis along 
the normal of the crystal plane and another axis along one direction of 
slip on that plane. With a, 8 referring to the subscripts of y's, the 
plastic shearing strain due to sliding on a plane with normal y, along the 
direction yg 


” 3 
Y = £€ a8 (5) 


where ug is the displacement along y,-direction. ‘This sliding system is 
here called af-slip system. With /,; and Js; denoting the cosines of 
the angles between y, and x,; ys and x,, respectively, the incremental 
strain components in x-coordinates due to de’’as, 


= (lial ig + sl jg) as. (4) 


Referring to strain coordinates defined as 2@:2 


For 1 = like de" ;; = 2 3) de’ 


For # 7 like = 2(Lialog + ig)de’ ap. 


There are as many de”’,s's as the active slip systems. de’’;; is caused by 
all the de’’,s's in the crystal. The resolved shear stress in the af-slip 
system due to stress tensor r,; to be critical 74s, may be expressed as (11) 


Tag = lad (9) 


where the repetition of subscripts « and 7 denotes summation from 
1 to 3. Writing this out in long form 


+ lad 22 + as + (lails2 + 12 
+ (lailas + + + 722. 
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This is a plane equation in stress space with 7,,;'s as coordinates and is 
called a yield plane. The normal to this yield plane has direction 
cosines with 7,,-coordinate proportional to the 7,; coefficient in Eq. 5. 
It is seen from Eq. 4, de’’,, is proportional to the same coefficient and 
hence normal to the yield plane. This is essentially the same as given 
by Bishop (12) showing single crystal yield surface consisting of yield 
planes which coincide with the plastic potential. 

With the consideration of elastic strain (13), the incremental stress- 
strain relationship of single crystals is uniquely determined. 


Associated Flow Rule of Polycrystals of Infinite Elastic Modulus 


A piece of metal consists of numerous crystals of various orientations. 
Consider a piece of unit cube subject to a uniform surface stress S,,, 
which is referred to as the macroscopic stress. The macroscopic strain 
increment is defined as 


dE;; = dau, + 1,du,)dS (7) 


where S is the area on the cube, /, is the cosine of the normal to dS with 
x-axis. The work done by the surface force S,; must equal the energy 
stored and that dissipated in the unit cube 


V. (9) 

kor infinite elastic modulus d/’,; and de’,; approach zero, the above 
gives 

S, dE” (10) 

For given dk”’,;, S,,ison the yield surface. Let S*,, be any stress vector 


within the yield surface, hence giving resolved shear stresses in all 
crystals in the cube less than the critical. 


| 
3 
: 
S, dk (r, de,dV. (8) 
. 
is Let single prime denote the clastic part ra 
& 
S* dak de’ ,;)d\ (11) 
(S;; — = f — > 0. (12) 
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This shows the validity of the associate flow rule for polycrystals of 
infinite elastic modulus. Bishop and Hill (10) derived the same equa- 
tion with S,; denoting £-7,,dA where dA is integrated over a unit area 
of any plane section. In their derivation no correlation between any 
component of the microscopic stress and any component of displace- 
ment over any plane section of unit area is assumed. With the present 
definition of S,; as the uniform stress over the faces of the unit cube, 
this assumption is not needed. 

For actual polycrystals with finite elastic modulus, S,,d/’;; may not 
equal /'7r,,de’;dV and Eqs. 10 to 12 do not follow from Eq. 9. Hence 
this proof is not applicable to actual polycrystals with finite elastic 
modulus. 


Associated Flow Rule for Polycrystals of Finite Elastic Modulus 


Actual crystals have finite elastic constants. For an aggregate of 
elastically isotropic crystals, any change of stress within the yield sur- 
face produces only elastic strain increment which is isotropic and homo- 
geneous and same for all crystals. The orientations of crystals in an 
aggregate are taken to be uniformly distributed in all possible planes 
and directions. The slip systems of all crystals cover all possible direc- 
tions on all planes. The initial critical shear stresses are taken to be 
the same for all crystals, hence a loading in stress space giving critical 
shear stress on any plane along any direction gives a point on the initial 
yield surface. Therefore the initial yield surface is a surface of maxi- 
mum shear stress. This has been shown by Lin (14) as 

> 

Js (4K? + J’,)(K? + 
27 

where J’, and J’; are the deviated stress invariants and K is a constant. 
When the stress vector reaches the initial yield surface, one slip system 
of the most favorably orientated crystal reaches critical shear stress and 
is ready to slide. Since the stress is uniform before sliding, the tangent 
plane to the yield surface of the aggregate at this point coincides with 
the yield plane of this slip system of this crystal. Infinitesimal sliding 
in this crystal occurs while other crystals remain elastic. Macroscopic 
plastic strain increment is due to the plastic strain increment of this 
crystal alone. This plastic strain increment of this crystal is perpen- 
dicular to the yield plane, hence di’, ; of the aggregate is normal to its 
initial yield surface. This shows that the initial vield surface coincides 
with the plastic potential. 

After sliding occurs in this most favorably orientated crystal, the 
stresses in the aggregate are affected. With body force F, per unit 
volume, the equilibrium condition (11) gives 


+ =0 (13) 
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where the subscript after comma denotes differentiation. The elastic 
strain equals total strain minus plastic strain, the plastic volumetric 
strain is taken to be zero. 


Nb, + (14) 


where 6,; is the Kronecker delta. and uw are Lame’s constants and @ is 
the volumetric strain. Substituting Eq. 14 into Eq. 13, 


= 0. (15) 


It is seen that — 2ye’’,, , is equivalent to F, in the strain field equation. 
Hence the strain field produced by slips in crystals is same as that caused 
by the equivalent body forces. From the strain field, the stress field is 
obtained. In a fine grained metal, the stress field due to e”;, ; in an 
interior crystal can be approximated by replacing e”’,,;,; by Fy = —2ue’’,;,; 
in an infinite medium. The stresses due to an applied force F, in an 
infinite medium have been shown by Kelvin (11) to be inversely propor- 
tional to the square of the distance from the point of the application of 
the force, the stress due to slip dies off rapidly with the distance from 
the slipped crystal and the main part of the aggregate is unaffected. 
Most crystals will still yield when the applied aggregate stress gives 
critical shear stress in one slip system of any of these crystals. Hence 
the yield surface is deformed in the vicinity where loading proceeds. 
Local expansion of yield locus in the initial stage of loading beyond 
elastic limit has been observed by Marin and Hu (15) and others. 
As loading proceeds, more crystals reach critical shear stress tag = T-. 
To keep these crystals from further slip, r.3 must be kept not to exceed 
r.. Hence dr,3 must belessthan zero. = 0 for each of these crystals 
gives a boundary plane tothe yield surface. These planes pass through 
the existing stress point and form a corner in the yield surface. In a 
combined tension and torsion test with S,;, denoting the applied aggre- 
gate tensile stress and S,, the aggregate torsional stress, the specimen 
is stretched beyond the elastic limit. Within the existing yield surface, 
all crystals are subject to the stress increments applied to the aggregate 
dr,; = dS,,;, although the total stress components may vary from one 
crystal to another due to slips. With a8 denoting a slip system reaching 
the critical value in some crystal, the boundary planes of the vield sur- 


face are given by 


S11 + + 12 = (), (16) 


This gives 


dS,» 
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This yield locus is similar to that predicted by the slip theory of Batdorf 
and Budiansky (16). However in the present paper, the stresses due 
to slips in the crystals are included so r,; # S,,;, while the slip theory 
assumes all crystals being subject to the same stress with 7,;; = Sj). 

Among all the af active slip systems of all the sliding crystals under 
tension, one gives maximum ratio of dS,,/d.S,. and one gives minimum. 
These two values of d.S;,/dS,. give the yield locus near the stress point. 
As stretching increases, more crystals slide, the range of a@’s giving 
Tas = Te increases and hence the range of d.S;,/dS,. given by Eq. 17 
increases, the maximum value increases and the minimum value de- 
creases. This indicates that theoretically as tension proceeds, the 
corner in the yield locus becomes sharper. This agrees well with the 
test data obtained by Naghdi, Essenberg and Koff (17). Tests by 
Budiansky, Dow, Peters and Shepherd (18) on thin-walled aluminum 
alloy cylinder in axial compression and then with additional combined 
compression and torsion show that the vield locus has a rather high 
curvature or sharp corner at the start of twisting. Here the yield 
surface is defined as the limiting surface in stress space, beyond which, 
incremental plastic strain, no matter how small, will occur. However 
in experiments, the yield locus is only observed when the incremental 
plastic strain reaches measurable magnitude. With improving ac- 
curacy in strain measurement, the corner in yield locus should be more 
apparent. 

After sliding starts in the aggregate, the stresses in the individual 
crystals consist of the uniform aggregate stress S,; and the non-uniform 
stress due to plastic slip denoted by 7,;,.. Let the resolved shear stress 
in aB-slip system corresponding to the aggregate stress S,; be 7,;, and 
that due to slips be 7.s,. The initial critical stress of a crystal is de- 
noted by 7., and the increase of critical shear stress due to hardening 
by Ar, 


In the active slip systems, 


Tap = ladgjSij + Tat, = Tie + Ate. (18) 


After the stress vector goes beyond the initial yield surface, a new 
yield surface is formed. Within this new yield surface the incremental 
stress produces no plastic incremental strain and the stress increment is 
the same for all crystals. Aras = ladg;AS,;. The limiting stress sur- 
face in the coordinates of S;; giving no plastic strain increment in all 
crystals, is the yield surface. The stress giving critical shear stress 
to a slip system in a crystal gives a boundary plane to the yield surface. 
Due to those crystals which have not slid previously the boundary 
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plane to the yield surface is 


Tos = bada)(Si; + QSi;) + Tea, = Tea (19) 


where the negative sign denotes sliding in the opposite sense. It is seen 
that the initial yield planes of these crystals are shifted by the amount 
Tag, [his shift is similar to kinematic hardening as proposed by Prager 
(19). Due to any of those crystals slid previously, the boundary plane 
to the yield surface is given as 


Tap = ladaj(Si; + ASij) + Tes, = (Teo + At-). (20) 


If the hardening is the same in different slip systems of the crystal, the 
initial yield polyhedron of the crystal is expanded by Ar, and shifted by 
the amount 7,3, in the stress space. This is a combination of isotropic 
and kinematic hardening as suggested by Hodge (20). Hardening 
Ar, may be different in different slip systems. If so, the different planes 
of the yield polyhedron move outward by different amounts by Ar,. 

The initial yield surface of the aggregate is bounded by the initial 
vield polyhedrons of all the individual crystals. After slips occur part 
of these yield polyhedrons of the individual crystals are translated, part 
are translated and expanded as discussed. The new yield surface is 
bounded by these new yield polyhedrons. ‘The initial yield surface of 
the aggregate does not shift nor expand as a whole, but the yield poly- 
hedrons of the individual crystals in the aggregate shift and expand 
differently. This corresponds to linear piecewise plasticity which has 
been studied by a number of distinguished investigators (21, 22). 

When the stress vector reaches a new yield surface, part of the 
crystals have the resolved shear stress equal to their critical values and 
are ready to start slide or continue slide. The remaining crystals of the 
aggregate have resolved shear stress less than the critical. An infini- 
tesimal increment of applied aggregate stress may cause some of the 
former crystals to slide by infinitesimal amount but does not give any 
finite stress increment so the latter crystals remain to have resolved 
shear stresses less than the critical and hence do not slide. The infini- 
tesimal slip in each of the sliding crystals gives incremental plastic strain 
vector normal to the yield plane of the corresponding slip system. The 
sum of the infinitesimal slips in those sliding crystals gives an incre- 
mental plastic strain vector lying within the cone of normals of the yield 
planes of these sliding crystals. Hence for the subsequent yield surface 
plastic potential exists and coincides with the yield surface. The rela- 
tive amounts of the increments of the resolved shear stresses in different 
sliding crystals depend on the direction of the incremental stress vector. 
Hence the direction of the incremental plastic strain vector dE”’,; due 
to the sum of incremental slips of all the sliding crystals depends on the 
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ratios of the stress components. Possible dependence of incremental 
plastic strain direction on incremental stress at a corner in the yield 
surface was first pointed out by Drucker (3). 

The above derivation is based on finite isotropic elastic constants 
which seems to be true for tungsten (23). Most metals have aniso- 
tropic elastic constants. The variation of elastic constants with direc- 
tions is rather large for a-iron, small for aluminum (23). With the 
consideration of anisotropic elasticity of individual crystals, the stresses 
vary from point to point in the aggregate even within the elastic range. 
The resolved shear stress in a crystal 7.3 is no longer equal to /,,3;7,;, 
but the sliding in this af-slip-system still gives an dE”’,; normal to the 
plane /,,/3;6;; = ladsjt:; = t-. Hence the plastic strain increment is not 
strictly normal to the yield surface. 

Bishop and Hill (10) have shown that the associated flow rule of 
plasticity applies to polycrystals with infinite elastic modulus. This 
paper shows analytically the validity of this flow rule for polycrystals 
of finite isotropic elastic modulus. However plastic potential deviates 
from the yield surface for aggregates of crystals of anisotropic elasticity. 
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ON THE PURE BENDING OF AN ELASTIC, PERFECTLY PLASTIC 
CURVED BAR IN THE STATE OF PLANE STRAIN* 


BY 


S. A. MURCH! 


SUMMARY 


The existing nonhardening, elastic-plastic solution to the problem of the pure 
bending of a curved bar is found to be valid only for a limited range of geometrical 
parameters, when the material is assumed to be compressible. In the complete treat- 
ment of the problem presented here, regimes of Tresca’s yield surface are encountered 
in which oe, is no longer the intermediate principal stress, rendering the problem 
statically indeterminate. Depending on the parameters of the problem, four types 
of solutions (with mutually exclusive domains of validity) are possible in which as 
many as six plastic zones with different stress fields may be present. Appropriate 
stress distributions are derived and complete results presented for a typical stage of 
deformation for one type of problem. In the limiting case of a completely plastic bar, 
it is found that the plane components of stress, the moment and the limiting position 
of the elastic-plastic boundaries are exactly those predicted by the rigid plastic solution. 


1. INTRODUCTION 


A solution to the problem of pure bending of a curved bar of non- 
hardening plastic material in the state of plane strain and employing 
the Prandtl-Reuss stress strain relations has been given by Shaffer 
and House (1, 2),2 where through the assumption of incompressibility 
the Mises yield condition takes on a special form which renders the 
problem statically determinate. Since for an incompressible material 
Mises’ and Tresca’s yield conditions (when the axial stress ¢, is the 
intermediate principal stress) are identical, the stress field, when Tresca’s 
yield condition is used, is also statically determinate. Furthermore, 
since the functional form of the displacement for an incompressible 
material may be found without recourse to the stress-strain relations for 
this problem, as shown in (3), it follows that the complete solution for 
an incompressible material is the same regardless of whether Mises’ 
or Tresca’s yield condition is used. 

Because the assumption of incompressibility does not alter the 
structure of the Tresca yield condition, the statically determinate stress 
field, given in (1), will still be valid in part for a compressible material,’ 
provided that ¢, is the intermediate principal stress in the plastic regions 


* The results presented in this paper were obtained in the course of research sponsored by 
the Office of Ordnance Research, U. S. Army, under Contract DA-04-200-ORD 171, with the 
University of California at Berkeley. 

! Assistant Research Engineer, University of California, Berkeley, Calif. 

2 The boldface numbers in parentheses refer to the references appended to this paper. 
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and remains so throughout the entire range of plastic deformation. \s 
will be shown presently, this requirement that ¢, be the intermediate 
principal stress places definite restrictions (depending on the value of 
Poisson's ratio and the geometrical parameters of the problem), on the 
range of validity of the solution given by Shaffer and House (1). Hence 
this solution, employing Tresca’s yield condition, if interpreted to hold 
for compressible materials, must be regarded as incomplete. 

The elastic-plastic problem as described in (1) begins with yielding 
at the inner radius r = a, forming a plastic region I (a = r = p,) with 
an outward moving elastic-plastic boundary p;. As deformation con- 
tinues, a second plastic region II (p, =r = 6) initiates at the outer 
radius r = 6, with an elastic-plastic boundary p, that moves inward 
until it meets p,;, as the curved bar becomes fully plastic. In Section 2 
of this investigation it is shown, however, that for any choice of the 
ratio 6/a of the outer radius to the inner radius of the bar, there exists 
physically admissible values of Poisson’s ratio v beyond which the state 
of stress in the plastic domain reaches a singular regime of the Tresca 
vield hexagon before the bar becomes fully plastic. Similarly, for each 
positive value of 0 < v < } there exists a value of the ratio b/a above 
which a singular regime of the yield hexagon is also reached in the range 
of contained plastic deformation. Hence, it is necessary to re-examine 
the character of the solution for the entire range of 0 < v = 3 and for 
all ratios b/a. Indeed, it is found that, depending on the values of b/a 
and », four different types of solutions are possible of which the statically 
determinate type given by Shaffer and House (1) is only one. This 
solution, however, constitutes (until the singular regime is encountered) 
the early stages of plastic deformation for all values of vy and 6/a. In 
later stages, the stress paths in the elastic domain enter a regular regime 
where ¢, is not the intermediate principal stress, thus rendering the 
stress field statically indeterminate. In Section 3, a stress field is 
derived for this regular regime and in Section 4 the solution, employing 
this statically indeterminate stress field, is found for one stage of the 
deformation and the solution for the fully plastic section is obtained as 
a limiting case. 


2. THE FOUR TYPES OF SOLUTIONS AND THEIR REGIONS OF VALIDITY 


In the state of plane strain and with reference to cylindrical polar 
coordinates (r, @, 2), we recall from the results given in (1) that the 
a, when the stresses satisfy the vield condition 


initial vield occurs at r 


oe — o, = 2k. (1) 

* Although in the determination of the stress field the assumption of incompressibility has 
not been introduced explicitly by Shaffer and House (1), a detailed examination of their work 
reveals this assumption is necessary if the solution presented is to be valid for all ratios of the 
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The nonvanishing components of the stress field in the ensuing plastic 
domain (zone |; a = r = p,) are given by 


2k 


r 
= kc (2a) 
a 


r 
= 1 + log-, (2b) 
a 


r 
(1 + 2tog" ), (2c) 
a 


where & is the yield limit in simple shear. The stress distribution given 
by Eqs. 2 remains a valid solution in zone I, only as long as a, is the 


Fic. 1. Tresca’s yield condition showing the stress profiles schematically. The dashed 
curve represents the state of stress at initial yield and the portions of the yield surface labelled 
\’B, B, BB’, DD’, D and DC’ correspond to states of stress in zones I, IIla, I11+, 11, [Va and 
IVh, respectively. Zone IV occurs in the tvpe D solution only 


intermediate principal stress. An examination of (2) in principal stress 
space (Fig. 1), however, reveals that as the elastic-plastic boundary 
p; moves away from r = a, the state of stress at r = p; approaches the 
edge (singular regime B) of the Tresca yield surface, described by 

o, — o, = 2k, (3a) 

regime Bs 
(3b) 
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The position p; = p; for which this singular regime is actually reached is 
found by substituting Eqs. 24 and 26 (which already satisfy Eq. 3a) 
into Eq. 3b; hence 


Pi v 
— = exp (4 
a | 1 — 2p ) 


The position r = pi, however, is not always attained by p,, since the 
values of 6/a and v may be such that the inward moving elastic-plastic 
boundary p, is met first. In this event the solution in (1) is valid, and 
the limiting position p, = p, is given by 


ab. (9) 
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Fic. 2. Regions of validity in the 6/a versus v plane for each of the four problem types 
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Fic. 3. Curved bar in stage III of plastic deformation (corresponding to problem types B, C, 


or D), showing schematically the relative positions of the elastic and plastic zones 


Thus, the inequality p; > p represents the condition that the stress field 
(Eqs. 2) be valid for a compressible material. By Eqs. 4 and 5, this 
becomes 


> exp 


(6) 
a 1 — 2p 


A plot of (6), when the equality holds, is seen in Fig. 2 (curve OA) to 
divide the plane of b/a versus v into two regions. For values of b/a and” 
v below OA the solution of Shaffer and House is valid; this will be re- 
ferred to, hereafter, as the type A solution. For values of b/a and » 
above OA, the state of stress at r = p,, reaches a singular regime and 
further expansion of the zone must be described by a new stress field, 
corresponding to the regimes of Tresca’s yield function encountered. 
The plastic zone which forms adjacent to zone I will be referred to as 
zone III (see Fig. 3). 

In a similar manner, from an examination of the statically determi- 
nate stress field in zone II (p, = r = 5) 


b 
= log -, 
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log 1, (0) 
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lying on the surface 


o, — a = 2k, (8) 


the state of stress at the elastic-plastic boundary p,, when considered 
in stress space, may be seen to approach the singular regime (Fig. 1) 
described by 

og = 2k, 


D: 2 ( 
regime D: ve (9) 


The reginie D is reached when 
= exp (10) 


but again depending upon the choice of b/a and v, p may be reached by 
p. before p; is reached in which case the stress field (Eq. 7) would re- 
main valid throughout the problem. It will be shown later (Section 4) 
that Eq. 5 is valid even in the presence of a plastic zone III, so that by 
Ikqs. 10 and 5 the criterion that the singular regime D, described by 
Kq. 9, never be reached is given by 


bh 


a 


2(1 
(11) 
exp 


Equation 11, when the equality holds, is plotted (curve CD) in Fig. 2; 
in the region above CD the stress at p, in zone II reaches the regime D 
of Tresca’s yield hexagon, and subsequent deformation must therefore 
be accompanied by the formation of a plastic zone IV (p, =r = ps) 
adjacent to zone IT. 

In order to classify the various types of solutions that may occur 
following initial yield and the formation of zone I, it is necessary to 
establish which of the zones II or III will appear after zone I. In the 
absence of zones III and IV the position of p,, when initial yield occurs 
at r = 6, may be obtained by placing p, = 6 in the relation between 


p, and p, 
2 b 
() tog ) = 1 + log (12) 
Pi pr ab 


Substitution of p; from Eq. 4 into Eq. 12 then gives the value of b/a for 
which zones II and III appear simultaneously. A plot of the resulting 
curve is shown by OB in Fig. 2; above OB zone III precedes zone II, 
while below OB and above OA zone II precedes zone III. Further 
analysis reveals zones II and III always precede zone IV so that only 
four types of solutions are possible; these are summarized in Table | 
and are to be identified with the corréspondingly labelled areas in Fig. 2. 
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Each of the four types of solutions mentioned above may be divided 
into a number of stages, corresponding to the number of plastic zones 
occurring ; that is, with the appearance of each plastic zone a new stage 
of the problem is entered. Since zone I is always the first plastic region 
to appear, the first stages of deformation for each of the four solution 
types are identical and the stress field is thus given by stage I of the 
solution given in (1), that is, the type A solution. Similarly, since zone 
I] is the second plastic region to appear in the type B solution, the stress 
field is given by stage II of the type A solution. For illustrative pur- 
poses the third stage of B, which is mathematically equivalent to the 
third stages of C and D is dealt with in the next section. Thesolutions 
in stages IT and IV of problem type D and in stage II of problem type © 
are not discussed, but each may be obtained in a manner similar to 
that of Section 4. 


I. 


Order of Appearance of Plastic Zone 


TABLI 


jx ot Solutio: 


A II 

B I I] 
| 
1) 


3. THE STRESS FIELD IN ZONE Ill 


Having established the ranges of b/a and vy for which the distribu- 
tions (Eqs. 2 and 7) are no longer capable of completely describing the 
plastic state of stress in all stages of deformation, we now proceed to 
derive stress fields that are valid in the plastic region II], corresponding 
to the new regimes encountered by the stress paths on the Tresca yield 
surface. Through the use of such stress fields, the complete solution 
for stage II of both the C and D problem types as well as for stage III 
of the B, C, and D problem types may be obtained. An examination 
of the loading paths (Fig. 1) in stage I of the type A solution (which is 
common to all solution types) shows that once the stress at r = p, in 
zone I has reached the singular regime B, further yielding may occur 
according to 


Og = 2k. 13) 


The presence of ¢, in Eq. 13 renders the stress field statically indetermi- 
nate and use of the flow rule 


Of of | 


| Oa, Oa 


(14) 


must be made in order to obtain the form of the stress field. 
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Since in the state of plane strain the total axial strain is zero, use 
of the flow rule (Eq. 14), together with Eq. 13, followed by an integra- 
tion with respect to time, gives 


(15) 


where prime and double prime refer to the elastic and plastic com- 
ponents of strain, respectively. Because the plastic strains are ex- 
pressible in terms of the elastic component of the axial strain, the total 
strains may be obtained in terms of the stresses through the generalized 


Hooke’s law and are given by 


1 
[o, — + a,) |, 


1 
+ O, te 20, + |, 


(). 


Substituting in Eqs. 16 for ¢, and a» from the yield condition (Eq. 13 
and the stress differential equation of equilibrium 


— (ra,), 17) 
or 


respectively, the appropriate equation of compatibility then becomes 


Oo, 


2(1 — + 6(1 — v)r + (1 — = 2(1 — »)RP, (18) 
or? or 
where P is an arbitrary function of time. The solution of Eq. 18 is 


+ (19a) 


where D, and D, are arbitrary functions of time, and 


Se 


€, 
( 1 6) as 
% 
= & 
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By Eq. 17 and the yield condition (13) 


D 
+ (1 — 53) + P, (19¢) 
2k 


(19d) 


It is found, however, that the stress field in zone III cannot be de- 
scribed completely by Eqs. 19 alone and still satisfy all of the con- 
tinuity requirements at the boundaries p, and p,. These conditions 
may be satisfied however, if zone III is divided into subregions Ila 
(op, =r = pz) and IIIb (pp =r = p,). If the state of stress in the sub- 
region IIIa is taken to lie at the singular regime B of the yield surface, 
described by Eqs. 3, then by Eq. 17, 


r 
= = log 5 (20a 
r 
= y (?0 
rb 1 + log a’ 20h) 


the stress field in the subregion III+ is given by Eqs. 19, while in the 


elastic domain the stresses are 


Po 
— 2 log +20, (21a) 
2 
(=) A+B (3 — 2 log ) + 2¢, (21d) 
2k r r 
| Po 


where A, B, and C are arbitrary functions of time.‘ 

It may be noted here that the stress field in plastic zone IV of stage 
IV of problem type D (which like zone III also consists of two sub- 
regions) is described by states of stress lying in the regular regime ED 
and the singular regime D of the Tresca yield hexagon (Fig. 1). These 
stresses may be derived in a manner similar to that used in obtaining 
Eqs. 19 and 20. 


‘See, for example, Timoshenko and Goodier (4). 
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4. THE STRESS FIELD IN THE PRESENCE OF ZONES I, I, AND Ill 


Since the solution for stages I and II of the problem type A corre- 
sponds identically to the first two stages of problem type B, we are 
motivated to complete the type B solution by considering stage III of 
this problem in some detail. In addition, with reference to Table 1, it 
may be seen that the presence of zone III together with zones I and II 
also typifies stage III of the problem types C and D, as well as the prob- 
lem type B. The form of the state of stress in the elastic and plastic 
zones, shown schematically in Fig. 3, may be summarized as follows: 
in zone I (a = r = p,) the stress field is given by Eq. 2 and satisfies the 
yield condition (Eq. 1); in zone IIla (p: = r = pz) by Eqs. 20, satis- 
fying Eqs. 3; in zone IIIb (p2 = r = p,) by Eqs. 19 satisfying Eqs. 13; 
in the elastic zone (p, = r = p.) by Eqs. 21; and in zone II (p, = r = 5) 
by Eqs. 7 which satisfy Eq. 8. 

From an examination of the stress field it is seen that the seven 
independent conditions that ensure the continuity of the stresses at the 
boundaries p,, p., and p,; together with the end condition 


are not sufficient to determine all of the nine quantities A, B, C, D,, 
D., P, M, p. and p,; for a given value of p; (which may be selected here 
as the time variable, since it increases monotonically throughout the 
deformation). Clearly, the stress field cannot be determined com- 
pletely without recourse to the displacement field which demonstrates 
further the statically indeterminate nature of the B, C, and D type 
solutions. It will be shown in Section 5 that six additional unknown 
functions of time are introduced in the displacement field along with 
seven independent continuity conditions at the boundaries p., p,, p1, and 
p>. Thus, considering the stress and displacement fields together, there 
are fifteen equations (nonlinear in p,, p;, and p;) in fifteen unknowns. 
The stress field, however, may be found to within the arbitrary function 
P and conclusions may be reached regarding the fully plastic solution 
without reference to the displacements; therefore a detailed treatment 
of the displacement field is left to the next section. 

The coefficients (A, B, C), appearing in the solution for the elastic 
zone, may be completely determined in terms of the positions of the 
elastic-plastic boundaries (p, and p,;) through the continuity of o, and 
oe at r = p, and by the satisfaction of the yield condition (Eq. 13) at 
r = p,, giving 


; 
‘ 
310 [J. F. 1. 
: 
(99 
ZZ) 
A 
; 


Oct., 1960. Pure BENDING or CurRvED Bar 311 


A= ‘| - 4(2 — v) + 4(1 — 2») log : | 
N pi 
1 b 
B= 6 + 4v — 2x? + 4(1 — 2r) log — | 
4 Po 
C= — 4v + x? + logx[4(1 — v) — 2x? —4(1 — 2v)logx] (23) 
b 
+ log — [2x? — 6 + 8» + 4(1 — 2r) log x}, 
where 


1+ 2x7 + 2(1 — 27) log x], 


and 
xX = po/ pi. 


Since ¢, is the intermediate principal stress in zone II, its continuity 
at r = p, is ensured by requiring the continuity of ¢, and o», since for 
the regime ED (Fig. 1) of the Tresca yield surface ¢, = v(a, + a). 
At r = p,, however, since a, is the intermediate principal stress in zone 
I11d, the continuity of ¢, is not implied by that of the plane components, 
so that there exists three independent continuity requirements on the 
stresses g,, a, and ¢,._ In this connection, it should be mentioned that 
equilibrium and the yield conditions only require that ¢, and the differ- 
ence — be continuous atr = p,._ If, however, and individually 
are allowed to be discontinuous, it is found through the compressibility 
equation that the displacement must also be discontinuous which is 
physically unacceptable. Hence, in addition to the requirement that 
(os — o.),, = 2k, employed above, we also require continuity in ¢, and 
a, at r = p,, so that by Eqs. 19 and 21 


{— sixtA + B(2 + — 2h log x) +.5,(2C — P)}, 


S 


(24) 


*2 
D, = —*— {sx*A — B(2 + 5, — 2s, logx) — 5,(2C — P)}. 


= 33 


At r = p2, since both of the stress fields (Eqs. 19 and 20) satisfy Eq. 13, 
there are only two independent stress continuity conditions, which may 
be imposed on and giving 


D, = 


D; 


: fi 
| 
S2 
P2 
1 +s: (log” P 
(2 >) 
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Equating the values of D; and D, obtained in Eqs. 24 and 25, leads to 
the following two equations describing the interaction of the boundaries 


Poy Piy Pr: 


— Ax’*s, + B(s, +2) — 2s,B log x — 2Cs, 


(s. log — 5,.P + 1), 


(26) 


— Ax’s, + B(s,; + 2) — 25,B log x — Ps; + 2Cs, 


= yt log + 1), 
( 


where y = p.,/pi. 

The above expressions correspond to the single relation (12) be- 
tween p, and p,; which holds on stage II of the type A solution. By 
setting p; = p: = p2 in Eqs. 26, where p; is given by Eq. 4, it is found 
that Eqs. 26 and 12 predict the same value for p,. This is to be ex- 
pected, since p; = p; = p2 represents the transition from stage II of a 
type B solution, where Eq. 12 is valid, to stage III, where Eqs. 26 hold. 
Allowing x to approach unity in Eqs. 26 leads to the result that 


lin y = 1, 
+] 


z 


so that in the limiting case of the fully plastic state, the plastic zone IIIb 
vanishes. This result is of particular significance since the stress field 
in this state is then given only by Eqs. 2, 7 and 20, so that the plane 
components of stress are seen to be identical to the plane components 
of stress determined by Shaffer and House (1), as the limiting case of 
an elastic-plastic problem (corresponding here to mode A) and identical 
also to the results of Hill (5), obtained directly using a rigid plastic 
material for the related problem of a bent sheet. Thus, by equilibrium 
the radial stress must be continuous in the limit so that by equating 
Eq. 7a with Eq. 20a, the limiting value of p,, p., and p,; is again given 
by Eq. 5. By Eq. 22 the bending moment per unit width of cross 
section, with the aid of Eqs. 24, 7b, 19c, 206, and 214, is found to be 


2M b\? 
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+ = [x? — 1 + 2(1 — 2p) log x] 
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x 1) (1 + s:tog® 
a 


— 1) (1 + s, log 


S; 


The limiting value of M in the fully plastic stage, that is, M = M** 
may be found from Eq. 27, following one application of L’ Hospital's 
rule, by allowing p;, p., 2 to approach p which leads to the result 


12 if, (28) 


| 


which is also in agreement with that of (1) and (5). 

Since, during Stage III of types B, C, and D solutions, both zones II 
and III are present (see table), the results (Eqs. 23, 24, 25, 26, and 27) 
of this section are valid for all three problem types; the expressions for 
the limiting case p./p; — 1, however, require separate treatment for 
problem type D due to the appearance of zone IV before the fully plastic 
state is reached. 


5. THE DISPLACEMENT FIELD 


For the sake of completeness and in order to determine the time 
function appearing in the solution for the stress field, the associated 
displacement field is determined in this section. In zones I, II, and IIa 
the displacements may be obtained through integration of the strain- 
displacement relations, where the strains are found through the equation 
of compressibility 


— 


E P= (o, + (29) 


and the appropriate condition of strain compatibility 


pie) 
or +2 or or ) 


where ¢,¢ has been taken equal to zero. Six arbitrary functions of time 
appear in the resulting expressions for u, and us. In the elastic zone 
and in plastic zone IIIa, since the total strains are expressible in terms 
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of the stresses by the generalized Hooke’s law and Eq. 13, respectively, 
no arbitrary time function other than that due to a rigid body motion 
term appears in the displacements which are found to be 


6 
(2 ) 2»)B log +K.—, (31a) 
kr r 
Gus ein 
= 4(1 — »)B6+ K, (31d) 


1 — 2p cos 


P+ + Ky (32a) 
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(32) 
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1+y 


1 — — 5.) v)(1 — ss) — 


(1 — s,)(1 4+ v) 1+ yp 


By the strain-displacement relations, Eq. 29 assumes the form 


1 — 2» 1 Aus 


r E (33) 


which is valid by the first of Eqs. 15 in zone IIla as well as the elastic 
zone. Since the form of us is such that the continuity of u» implies that 
of t it is seen from Eq. 33 that the continuity of u, at r = p; is a 
consequence of the continuity of the stress field and need not be further 
imposed. The six unknowns and P may thus be found from the re- 
maining seven continuity conditions in displacements at the boundaries 
pi, Px, p, and p,. Imposing these requirements, the radial component 
of displacement in zones I, II, and IIla becomes, respectively 
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and the tangential component of the displacement in all zones including 
the elastic zone and IIla is given by 


4(1 — »)Bo — (37) 


kr 


with the particular result 


1 — » 1 


The time function K = K, = K, is the same for each zone and may be 
evaluated by constraining a point on the 6 = 0 axis. This is most 
conveniently done by setting® u,(r, 6, t) equal to zero at r = p, @ = 0; 
that is, 

u-(p, 0, t) = 0, 


since r = p is the last remnant of the elastic zone. Thus K may be 
found from Eq. 31a in all stages of the problem 


6. DISCUSSION AND RESULTS 


Although the regions of validity of the A, B, C and D types solution 
are mutually exclusive (Fig. 2), it is of interest to consider the error 
introduced by regarding the type A solution as valid for all values of 
0 < v = 3 and the ratio b/a (in this connection see footnote 3 in the 
introduction), including the solution discussed in Section 4. For a 
given choice of b/a this is equivalent to assuming that Poisson's ratio 
lies to the right of curve OA in Fig. 2. It is found that for all practical 
purposes, the presence of zone III has little effect on the value of the 
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Fic. 4. Comparison of the value of po predicted by the type C solution (shown solid) 


for b/a = 20, v = 0.3 with the value predicted by the type A solution (shown dashed) for 
the same value of p,. 


bending moment even for moderately large values of b/a; for example, 
for b/a = 20, v = 0.3, the bending moment obtained from Eq. 27 
differs by less than one-tenth of one per cent from the value predicted 
by the type A solution which is actually only valid for » > 0.375. The 
predicted positions of the elastic-plastic boundaries however, differ 
considerably as may be seen from Fig. 4. 
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THE EVAPORATION OF THORIUM METAL 
BY 


DANIEL L. GOLDWATER AND WILLIAM E. DANFORTH 


Some years ago we wished to know the rate of evaporation of bulk 
thorium metal. C. Zwikker in 1928 included the results of his own 
unpublished measurements in a review article (1).'_ Values quoted by 
other authors (2, 3, 4, 5, 6) for heat of evaporation of thorium metal or 
its rate of evaporation, seem to refer back to this work. 

In the hope of quickly checking the published data we set out to 
measure the rate of evaporation of thorium from a tungsten dish heated 
in vacuuo by RF induction. Measurements were made by observing 
the weight loss of the specimen and also by finding the gain of weight of 
a nearby collector. Appropriate precautions were taken to define the 
evaporating area, obtain true temperature values, etc. In the earliest 
experiments, the thorium was exposed to total pressures somewhat 
above 10-§ mm. of Hg during the heating cycle. At temperatures 
around 1500° C., data from these experiments were in fair agreement 
with the previously published data. At lower temperatures, however, 
very erratic data were obtained, the specimen sometimes even showing 
a small gain of weight. Furthermore, it was noticed that the thorium 
had a golden patina when the specimen was removed from the vacuum 
system for weighing. 

In subsequent experiments, the thorium was maintained in a 
vacuum always better than 2 X 10-* mm. of Hg during the entire cycle, 
and care was taken that the specimen had completely cooled before the 
vacuum system was let down to air. The thorium was then found to 
keep a silvery lustre, and the evaporation data, although they still 
exhibited considerable scatter, were consistently 5 to 10 times higher 
than the previously published values. A comparison of earlier data 
and our data is shown in Fig. 1. 

It may be presumed that the patina and the low evaporation rates 
obtained at the ‘“‘higher’’ pressures were both due to surface oxidation 
of the thorium. Thorium oxide has a very low evaporation rate at the 
temperatures concerned. The slope of either set of data is entirely 
inconsistent with the heat of evaporation of thorium oxide (7). Thus 
it appears that a surface oxide forms at “higher’’ pressures which in- 
hibits, but does not entirely block the evaporation of the underlying 
metal. Even after improving the vacuum conditions, we occasionally 
observed small dark specks on the thorium surface which were probably 
residual islands of oxide. 


' The boldface numbers in parentheses refer to the references appended to this note. 
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The evaporation data may be represented by the form 
r = Aexp (—6b/T7), 
where the values obtained from our data are: 


A = 2.3 10° g-cm~?-sec™, b = 68,000° K. 


10°/T 


kG. 1. Specific evaporation rate of thorium vs. reciprocal temperature. Solid points 
loss of weight method. Open circles—collector method. Dotted line represents earlier 


published values (see 6). 


The heat of evaporation for thorium as derived from these data is about 
140 Kcal-mol-'. With some difficulty, measurements could be carried 
out on thorium specimens of greater purity and having still cleaner 
surfaces. In the absence of such measurements, the data presented 
here are believed to be an improvement on those previously available, 
and appropriate for many technological purposes. 
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Cosmic Rapio Waves, by I. S. Shklovsky, 
translated by Richard B. Rodman and 
Carlos M. Varsavsky. 444 pages, illustra- 
tions, 6 X 93} in. Cambridge, Harvard 
University Press, 1960. Price, $12.50. 


Since the first discovery of radio waves of 
extraterrestrial origin by Karl Guthe Jansky 
in December 1931, this new branch of science 
has made remarkable progress, particularly in 
the last ten years. It has opened new “win- 
dows” enabling us to “peep” through and ex- 
plore the Universe. The fascinating ‘views’ 
through these “windows” have been respon- 
sible for many discoveries of fundamental 
cosmological importance that constantly en- 
large our concepts about the Universe 

It has now been definitely established that 
radio waves are emitted from the Sun's 
gaseous atmosphere from our Milky Way 
Galaxy, from galaxies in collision at far away 
distances in space, and from numerous in- 
visible “heavenly"’ bodies. As the name im- 
plies, Dr. |. S. Shklovsky's book Cosmic Radio 
Waves deals with radiation only from objects 
far away from our solar system. This limi- 
tation has resulted in a detailed study of these 
waves. 

The book is divided into six chapters. The 
first chapter starts with the concepts and 
definitions of the basic quantities which are 
measured in radio-astronomical investiga- 
tions. Then, characteristics of receivers and 
antennas, the arrangements of antennas to 
form interferometers, and efforts to increas« 
the resolving power are described. It is 
stressed at the end of the chapter that the 
radio-astronomical methods form an essential 
supplement to optical methods, and real prog- 
ress depends upon the close cooperation be 
tween the two disciplines 

In the second chapter, surveys of the in- 
tensity distribution of cosmic radio waves at 
different frequencies made during most of the 
major investigations up to the present time 
are presented. Frequency dependence of the 
intensity of cosmic radio waves and the spatial 
distribution of the components of radio emis- 
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sion are discussed. From the analysis of the 
above data, the important conclusions—that 
sources of cosmic radio waves form a galactic 
disc with fairly strong concentration to the 
center, and other sources of galactic radio 
waves form a gigantic, almost spherical “halo” 
are reached. Numerical data of observa- 
tions on radio emission from discrete sources 
their angular extent, and optical identification, 
are next presented. It is pointed out that 
the term “radio star,’ often applied to dis- 
crete sources in general, is a misnomer since 
all discrete radio-wave sources identified up 
to the present time are associated with neb- 
ulae, either galactic or extra-galactic, and 
since not a single source has been identified 
with a star. Also it is shown that the spatial 
distribution of the stronger and weaker dis- 
crete sources of cosmic radio waves is opposite 
to that observed for the stars. The second 
chapter is full of wealth of the basic observa- 
tional material, presented in a semi-historical 
survey, and will serve as a useful ready refer- 
ence for the research workers in this field. 

Chapter 3 deals with problems of the phys- 
ical processes responsible for observed radio 
emission from the Galaxy. Arguments that 
the flat component of galactic radio waves 
is a consequence of emission by ionized in- 
terstellar gas are presented. Evidence for the 
interstellar magnetic fields comes from the 
polarization of light from some stars. Em- 
phasis is laid on the conclusion that the mag- 
netic-deceleration theory for the non-equi- 
librium radio emission of the spherical com- 
ponent can satisfactorily explain all the 
observational data. 

Monochromatic galactic radio waves are 
discussed in Chapter 4. Observations of the 
emission and absorption of the 21-cm. line 
of hydrogen put in the hands of investigators 
the most powerful tool for studying the spiral 
structure of the Galaxy, and already ex- 
tremely valuable information about the phys- 
ical properties of the interstellar medium, the 
dynamics and the structure of the Galaxy has 
been gathered. It is pointed out that the 
analysis of the 21-cm. line profiles provides 
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a new possibility for studying the random 
motions of clouds of interstellar gas. Pos- 
sibilities of finding other lines of deuterium, 
nitrogen and also molecular lines are outlined. 
This chapter gives a very good account of 
“radio astrospectroscopy.” 

The connection between radio-astronomy 
and the origin of cosmic rays is the subject 
of Chapter 5. Optical identification of the 
powerful discrete sources of radio waves with 
nebulae which are located at the sites of super 
novae explosions is very suggestive. In order 
to answer whether all galactic discrete sources 
of nonequilibrium radio emission be regarded 
as remnants of super novae explosions, the 
author carried out an investigation by seeking 
the information on such explosions in the 
ancient astronomical records of many coun- 
tries. He has tried to identify discrete radio 
sources at such sites, and reached the con- 
clusion that a super nova explodes every 30 
to 60 years, on the average. The frequency 
of such explosions is of extreme importance 
in questions of cosmic-ray origins. After 
looking into the nature of the radio emission 
by the ejecta of super nova, plasma-oscillation 
mechanism is discarded and the mechanism 
of deceleration radiation by relativistic elec- 
trons in interstellar magnetic fields to explain 
the peculiarities of the radiation from the 
ejecta of super nova is accepted. One full 
section is devoted to the optical continuum 
and polarization of the crab nebula, which 
lends support to the deceleration radiation 
mechanism. The last section of this chapter 
looks into the question of the origin of cosmic 
rays. Here the problem is to identify sources 
in which they arrive, and of acceleration 
mechanism in these Arguments 
against the proposd cosmic induction ac- 
celeration mechanisms are put forth. De- 
celeration radiation from the relativistic elec- 
trons shows the existences in the interstellar 
space of particles having energies from 10* 
to 10° ev, which is supposed to be the electron 
component of the primary cosmic rays. The 
spherical system formed by the sources of 
the general field of nonequilibrium galactic 
radio emission outlines the domain of the 
space within which primary cosmic rays are 
localized. It is regarded that super nova out- 
bursts maintain the dynamic equilibrium of 
the cosmic rays in our galaxy. Thus radio- 
astronomy has opened a completely new ap- 
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proach in our understanding of the nature 
and the origin of cosmic rays. 

The last chapter deals with extragalactic 
radio waves, and the nature of the radio 
galaxies. In the end the intimate relation 
between radio-astronomy and cosmology is 
stressed. 

This book outlines the progress in the field 
up to the 1958 Paris Symposium on Radio 
Astronomy. In addition the author has 
drawn attention to a few selected papers 
published since that time. In the appendix, 
these papers are briefly summarized and each 
one is keyed to the text of the book dealing 
with the same topic. This has made the book 
quite up to date. References are grouped 
under separate years covering the period 1846 
to 1960. This chronological order emphasizes 
the evolution of the field. In addition a 
general bibliography is provided offering a 
list of further reading. An Index is provided 
which includes both authors and subjects. 

This book is very useful to the active re- 
search workers in the field of radio-astronomy, 
because it focuses attention upon the major 
problems in urgent need of further study. 
Also it should be read by those who are in- 
terested in probing the mysteries of our 
universe. 

V. R. Pornis 
Bartol Research Foundation 


X-RAY ABSORPTION AND EMISSION IN ANA- 
LYTICAL CuHeMistTRY, by H. A. Liebhafsky, 
H. G. Pfeiffer, E. H. Winslow and P. D. 
Zemany. 357 pages, illustrations, 6 x 9} 
in. New York, John Wiley & Sons, Inc., 
1960. Price, $13.50. 


This book has been written by a team of 
authors who themselves have contributed ex- 
tensively to the recent developments in X-ray 
emission spectrography (X-ray fluorescence 
analysis) and X-ray absorptiometry. This 
is well demonstrated by the author index of 
810 entries of which 83 refer to the authors. 
The book has been written for the analytical 
chemist or any other worker whose primary 
interest lies in obtaining an analytical result. 
This frequently stated intention of the authors 
has to a great extent influenced the presenta- 
tion of the subject matter. Any newcomer 
to the field will find all the necessary informa- 
tion which will enable him to carry out a 


successful analysis. 
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In 11 Chapters and 7 Appendices, all well 
documented by some 450 references, there 
are treated the origin and properties of X-rays, 
X-ray absorptiometry with polychromatic and 
with monochromatic beams, X-ray emission 
spectrography and the determination of film 
thickness by X-rays. In the introductory 
chapters dealing with the physics of X-rays, 
the presentation follows closely the historical 
development. Here the aim of the book is 
most clearly visible; the physical relations 
are stated, if possible, in a form suitable for 
numerical evaluation without going into the 
details of their derivations. The analytical 
chemist to whom the book is addressed can 
either take these relations for granted, or he 
can read the original literature to which 
frequent reference is made. 

In the following chapters, which deal with 
the spectrochemical analysis with X-rays, an 
excellent account is given of the methods and 
tools used at the present time. By frequent 
reference to practical problems and_ their 
solutions, the scope and limitation of every 
analytical method mentioned are convincingly 
demonstrated. In order to obtain the same 
information without the help of this book, 
many, often obscure, journals would have to 
be consulted. Here, it seems, lies the im- 
portance of this publication. It collects in a 
well organized manner all the important as- 
pects of the chemical analysis by X-ray tech- 
niques excluding X-ray diffraction. 

Throughout the text the authors aim at 
systemizing the terminology which, not sur- 
prising for a rapidly developing field, is often 
confusing. However, the distinction between 
a spectrometer, which measures the wave- 
lengths and intensities in X-ray beam gener- 
ated by a tube, and a spectrograph, which 
measures the same quantities in a beam gen- 
erated by an unknown sample, appears rather 
artificial. The term intensity, which plays 
an important part throughout the book, has 
apparently resisted consistent interpretation 
The statements on page 6 do not give a clear 
definition of the term intensity, especially 
when considered in conjunction with the men- 
tioning on page 8 that the integrated intensity 
can be measured in watts. However, these 
are minor points which by no means distract 
from the extreme usefulness of this book. 

F. R. L. ScHONING 
The Franklin Institute Laboratories 


THe THEORY OF LINEAR ECONOMIC MODELS, 
by David Gale. 330 pages,6 * 9in. New 
York, McGraw-Hill Book Co., Inc., 1960. 
Price, $9.50. 


This book is an excellent example of what 
one hopes is a healthy trend in the application 
of mathematics to the social sciences. Not 
too many years ago the typical econometrician 
or psychometrician was an economist or 
psychologist who utilized in pretty much of 


a cookbook manner ‘“‘appropriate”’ statistical 
or mathematical techniques with little atten- 
tion to how they were developed and the 
assumptions which underlay them. There 
have been, of course, notable exceptions, as 
there usually are to such generalizations, and 
to them we owe the development of these 
fields as they are today. 

Now Professor Gale has written a book in 
mathematics which is of primary concern to 
economists, but which few students of eco- 
nomics will be able to read without consider- 
able effort. Mathematicians from other areas 
will not find it easy either. This situation 
arises because the modern mathematics which 
is basic to economic models is just now finding 
its way into the undergraduate curriculum, 
and because mathematics is enlarging its role 
in all the scientific disciplines faster than it 
is being made a part of most scientists’ 
training. While these remarks may imply a 
criticism of this volume from one point of 
view, this reviewer prefers to take the position 
that “growing pains’’ cannot be avoided in 
the development of good quantitatively ori- 
ented economists, and that the author's efforts 
toward this goal are admirable. He sets a 
standard even though it may be a while before 
there are many who will meet it. 

A brief outline of the contents serves to 
indicate the scope of the book. Chapter 1 
presents examples, definitions, and principal 
theorems of linear programming. Economic 
equilibrium which is used throughout as a 
unifying concept is also introduced here. 
Chapter 2 developes the needed real linear 
algebra, from vectors and matrices to convex 
sets and polytopes, for the remainder of the 
book. Chapters 3, 4 and 5 discuss, respec- 
tively, the theory of linear programming, the 
simplex method of solution with computa- 
tions, and the treatment of integral-valued 


problems such as the optical assignment of 


3 
4 
3 
3 
: 


Oct., 1960. | 


persons or machines. Chapter 6 is devoted 
to game theory and covers numerous ex- 
amples, definitions, and elements of the 
theory. Chapter 7 takes up the solution of 
matrix games showing their relationship to 
linear programming. Chapters 8 and 9 take 
up the more specialized economic models for 
exchange and production. The former dis- 
cusses the simple linear exchange model, dy- 
namic theory, and the place of price equilib- 
rium in exchange models. The latter chapter 
covers a simple and general linear model, and 
Neumann models. 

The exercises are numerous and varied in- 
cluding proofs, literal problems, and numerical 
examples. They are effective in promoting 
understanding at the levels of intended cover- 
age. The typography and format make for 
good legibility. Assumptions are spelled out 
so that their effect on the development is 
apparent. 

Two minor points which this reviewer, due 
to his statistical bias, would like to see in- 
cluded are some mention of stochastic input 
and its effect on the models, and adherence 
to what appears to be a helpful convention in 
books on multivariate analysis—the use of 
bold-faced type to denote vectors and 
matrices. 

It is hoped that the author will be suffici- 
ently enthused, as non-linear models become 
more thoroughly investigated, to write a com- 
panion volume on them which wiil be as good 
as the present text is in its chosen field. 

W. Jay MERRILL, JR. 
General Electric Co. 


CHEMICAL Periopicity, by R. T. Sanderson. 
330 pages, illustrations, 8 X 10} in. New 
York, Reinhold Publishing Corp., 1960. 
Price: $9.75 (college) ; $11.75 (trade). 


This is an unusual book, with an unusual 
approach to chemistry. In his introduction 
the author wonders ‘‘whether civilization in 
general might not ultimately gain more if 
scientists would affect a moratorium on all 
new experimentation for a year or two and 
spend time in leisurely contemplation of the 
facts already discovered.’ It seems that the 
author has done just about that, focussing 
his attention primarily on inorganic chemistry. 
And the conclusions and tables show that his 
contemplation was far from leisurely 
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Modern chemistry is (or ought to be) based 
on the periodic table of the elements and on 
the knowledge of the electronic structure of 
the atoms. The author presents chemistry 
under this light, and succeeds in showing the 
internal and logical consistency of the whole 
field of chemical reactions and affinities. Al- 
though our understanding of the chemical 
bond in molecules and in solids is far from 
established, it has progressed considerably 
since the advent of wave mechanics forty 
years ago. Enough progress has been made 
to explain not only why atoms combine but 
also to justify and sometimes even predict 
many of the properties of the combination 
(melting points and heat of fusion, boiling 
points, conductivity 
Cx. 

After an introduction to the periodic table 
(with very pertinent historical remarks) and 
to chemical combination, the author analyzes 


crystal structure, 


the principles of coordination chemistry and 
the properties of the individual elements. Al- 
most the whole second half of the book is 
devoted to the binary chemistry of non- 
The book more than 30 
of the properties of the 


metals. contains 
“periodic tables” 
elements that will prove invaluable. It is 
probably the first time that the properties of 
the elements have been presented in such a 
unified, extensive and logical fashion. The 
author is to be congratulated for this beautiful 
synthesis. 

It is regrettable that the quantum-mechan- 
of the is rather 


It need not have increased seriously 


ical basis chemical bond 
sketchy. 
the size of the book, since semi-quantita- 
tive arguments are possible even in wave 
mechanics. 

With this reservation expressed, this reader 
feels that this book will contribute enormously 
to improve the teaching and the understand- 
ing of chemistry, and to develop a ‘more 
reliable chemical intuition."’ It extends the 


concept of electronegativity and proposes 

several new interpretations. The large format 

of the book is well suited for the tables, and 

the publishers should be congratulated for the 
beautiful and easy to read printing. 

HENRI AMAR 

The Franklin Institute Laboratories 


and Temple University 
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ConFLict oF INTEREST & FEDERAL SERVICE, 
a report by the Special Committee on the 
Federal Conflict of Interest Laws of the 
Association of the Bar of the City of New 
York. 336 pages, 6 * 9} in. Cambridge, 
Harvard University Press, 1960. Price, 
$5.50. 


This book touches upon one of the most 
important problems which faces the Federal 
Government—‘‘Why is the Federal Govern- 
ment so inept as an employer that the most 
capable men and women will not undertake 
or retain Federal employment except in an 
emergency?” “Why does the Federal Gov- 
ernment have such a genius for creating an 
accumulation which 
many of the most able employees are not 


of minor annoyances 


willing to endure?’’ These questions have 
many answers. 

Certainly a part of the problem is that the 
technique for hiring and promoting personnel 
in Federal employment has been highly un- 
satisfactory. Thinking in terms of engineers 
and scientists, Federal employment and pro- 
motion has encountered great difficulty in 
recognizing the person who has learned ‘‘things 
the hard way,"’ but has extraordinary inge- 
nuity, who might be a leading research engineer 
in a corporate organization. Furthermore, 
Federal employment and promotion have 
been very naive about people with paper 
qualifications who lack real ability. 

If Federal employment has encountered 
serious problems in hiring and promotion, 
the difficulty in firing an incompetent Federal 
employee is even more serious. Those who 
have worked for the Government 
employees 


will re- 


member how incompetent are 
transferred from place to place. 

The genius for annoyance of Federal em- 
ployees is based upon an accumulation of little 
things. One practice is to give a worker a 
job to do without providing him with the 
necessary tools to do it—for example, without 
permitting him adequate secretarial help, or 
use of a long distance telephone, or requiring 
him to travel, and then making a series of 
unrealistic and inequitable rulings which de- 
prive him of proper reimbursement for travel 
expenses. 

The Association of the Bar of the City of 
New York is to be congratulated on their 
recognition of the great importance of this 


Book REVIEWS [J. 


problem and their painstaking, interesting and 
intelligent consideration of one phase of the 
problem. They correctly evaluate the present 
situation: 
The United States has been thrust into 
a position of world leadership encumbered 
with the staffing system, traditions, and 
attitudes toward public service of the 
nineteenth-century laissez faire govern- 
ment of a second class power. (Page 142) 


They recognize that the present morass of 
confusing statutes regarding conflict of in- 
terest is a serious obstacle toward employment 
of competent personnel by the Federal 
Government. 

They also recognize that this is a particu- 
larly important problem to the scientist. Ac- 
cording to the report: 

The old caricature of the scientist as a 

kind of combined Dr. Frankenstein and 

Mr. Chips has been replaced by a public 

image of the scientist as a hard working 

citizen of learning and probity, engaged 
in vital tasks and worthy of great respect. 

(Page 167) 


Actually the conflict of interest problem 
from the standpoint of scientific ethics should 
not be any different for the federally employed 
than it other 
scientist. Any man who is acting as a con- 
sultant, or who represents a company doing 
consulting work is urged to read particularly 
the illuminating comments regarding the ob- 


scientist would be for any 


ligation of the scientist in respect to conflict 
of interest. 

The Report realizes that it is completely 
unrealistic in most cases to impose on an ordi- 
nary man entering Federal employment the 
frightful price exacted of Charles Wilson when 
he became Secretary of Defense. You will 
remember that he was required to dispose 
completely of all of his General Motors stock. 

The conclusion reached by this Report is 
that we should make a completely fresh start 
in conflict of interest matters of the Federal 
Government. This is not the place to discuss 
the merits of the particular statutory changes 
which are proposed, although they cannot 
fail to be a great improvement over the 
existing statutes 

Contrary to the usual rule when a report 
“that it takes a 
this Report 
is beautifully written from start to finish 


is written by a committee 


committee to understand it” 


4 
; 
: 
: 


Oct., 1960.] 


A great deal of progress will be made if 
this book will bring home to those responsible 
for determining the of Federal 
employment, the fact that a fresh start is 
needed not only in the matter of conflict of 
interest, but on all phases which are concerned 
with hiring, firing and determining conditions 
of work of Federal employees. 

JoserH Gray JACKSON 
Wm. Steell Jackson & Sons 


conditions 


FUNDAMENTALS OF GaAs Dynamics, (VOLUME 
or HiGH SpEED AERODYNAMICS AND JET 
PROPULSION), edited by H. W. Emmons. 

Princeton, 

Price, 


749 pages, diagrams, 6 * 9 in 
Princeton University Press, 1958. 
$20.00. 


This book is Volume III of a 
volumes devoted to high speed aerodynamics 
and jet propulsion. The book opens with a 
Foreword by Von Karman acknowledging the 
support for publication of the series by the 


series of 


Armed Forces. 

The contents cover a comprehensive ana- 
lytical discussion, including basic experimental 
data on various aspects of gas dynamics. 
While the separate contributions on subject 
matter are listed in alphabetical order, by an 
outstanding group of scientists, they are 
uniquely edited to form a continuous and 
excellently coordinated and balanced discus- 
sion. Thus the work can be used not only 
as a reference work but likewise can be intro- 
duced as a basic graduate text as a research 
background for extending our knowledge on 
gas dynamics. 

The opening section, “A,” by H. S. Tsien, 
covers the general equations of gas dynamics. 

Section “B"’ by E. Crow is a one-dimen- 
sional treatment of gas dynamics. An excel- 
lent discussion of thermodynamic properties 
including compressibility effects, etc. and 
covering extensive data on air, precedes the 
discussion on flow relations. A critical dis- 
cussion of shock waves, extending to the flow 
in real nozzles, flow in diffusers, long ducts, 
is given 

Crow's contributions cover further exten- 
sive analysis on various types of flow, in 
ducts and nozzles, including the effects of heat 
transfer, specific heats and gas imperfections, 
of particular value in modern research in var- 
ious applications of gas dynamics. 
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Section “C"' by A. Kantrowitz is a very 
interesting discussion on non-steady flow, in- 
cluding the building up of the wave equation, 
with the effect of entropy changes. The anal- 
ysis extends to large amplitude isentropic 
waves. Stability of normal shock waves is 
followed by a discussion of very intense ex- 
plosions with a clear physical picture of such 
phenomena, including examples in rockets, 
and the “Kadenacy Effect"’ in internal com- 
bustion engines. 

Section “D" by Wallace D. Hayes dis- 
cusses the effect of discontinuities, with its 
relations to shock wave analysis. 

Section “E" by H. Polacheck and R. J. 
Seeger is on the important subject of shock 
wave interactions. The work includes a group 
of unusual photographs of actual effects of 
shock wave phenomena, such as the regular 
refraction, air-CO, and irregular refraction, 
air-CH, matters, with actual photographs of 
spinning detonation are shown. Along with 
this work the authors include an extensive 
analysis of refraction shocks. 

Section “F"’ by H. Guyford Stever studies 
the effects of condensation phenomena in high 
speed flows, and the general kinetics of con- 
densation, including experimental data. 

Section ““G' by Von Karman is a basic con- 
tribution to general aspects of combustion. 
In this, he reduces the complicated phenom- 
ena of combustion to basic parameters starting 
with the continuity equation, the momentum 
equation and the energy equations as con- 
Then, intro- 
ducing the concept of diffusion as the mass 
flux and a parameter for accounting the prog- 
ress of with chemical kinetics 
as a function of local temperature, pressure 
and chemical composition of mixtures, he 


stants for the entire process. 


combustion 


arrives at generalized equations for the basic 
phenomena of combustion. 

Section “G" is continued with the separate 
contribution by the editor, H. W. Emmons, 
on flow discontinuities associated with com- 
bustion. This is a scholarly discussion of the 
complicated subject of flame front extending 
the fundamental equations of heat transfer, 
with relating to the 


boundary conditions and the inadequacy of 


particular emphasis 


property reaction rates for the various mix- 
tures. However, the author covers a very 


illuminating and valuable discussion relating 
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to flame fronts, both oblique and general and 
steady flames. 

The work concludes with Section G-10 by 
G. 1. Taylor and R. S. Tankin. First the 
Chapman-Jonquet theory of detonation is 
discussed, followed by a thorough analysis of 
general theory of deflagration and detonation, 
including a very interesting discussion on 
spinning detonation in gaseous mixtures. 

RUPEN EKSERGIAN 
The Franklin Institute Laboratories 


THEORY OF INERTIAL GUIDANCE, by Connie 
L. McClure. 
in. New York, Prentice-Hall, Inc., 
Price, $12.00. 


340 pages, diagrams, 6 = 9 
1960 


There has been a need for some time for 
an introductory presentation of this popular 
area of research and especially for relating it 
to its classical dynamic This 
task is undertaken in this book by C. L 
McClure who has been an active investigator 


foundations. 


in this field for many years. 

The heart of the subject is dealt with in 
chapters seven to ten (about two-thirds of 
the book) 
lized platforms is first presented and followed 


The preliminary theory of stabi- 


by various single-channel descriptions of 

rhe presenta- 
The last 
chapter on Terrestrial Inertial Guidance Sys- 


their 


lerrestrial Inertial Systems. 


tion is lucid and authoritative 


tems and simulative evaluation is 
excellent. 

In the first six chapters, devoted to a review 
of the fundamentals of dynamics and gyro 
dynamics, it is regrettable that a rather 
clumsy notation has been used making it 
harder to follow otherwise simple arguments 
The appendices and bibliography will be 
found very useful. 

HENRI AMAR 
The Franklin Institute Laboratories 


and Temple University 


BIOCHEMISTRY OF STEROIDS, by Erich Heft- 
mann and Erich Mosettig. 231 pages, dia- 
grams,6 X9in. New York, Reinhold Pub- 
lishing Corp., 1960. Price: $5.75 (College) ; 
$6.90 (trade). 


The development of the biochemistry of 
steroids is more recent and has followed the 


structure elucidation of several classes of 


steroids by the organic chemists. The publi- 
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cation of “Chemistry of Natural Products 
Related to Phenanthrene"’ by the Fiesers in 
1936 and its two subsequent expansions, the 
latest resulting in “Steroids” in 1959, traces 
the rapid evolution of this field. The present 
work intends to supplement ‘“‘Steroids’’ and 
describes distribution, biosynthesis, metabo- 
structural 
specificity and analytical determination of 
each class of steroids. 


lism, physiological significance, 


The classes of steroids are arranged in the 
Ster- 
ols, steroid sapogenius and alkaloids, cardiac 
glycosides, bile acids and sex and adrenal hor- 
The subject 
matter is presented in textbook style and in 
All references (774) are 
tabulated at the end of the book and are very 


sequence of their biogenetic relations. 


mones are treated in this order. 
an interesting way. 


conveniently classified according to classes of 
steroids. The illustrations are very clear. 
Nomenclature and conventions in the repre- 
sentation of structure formulae correspend to 
those used in “Steroids.” 

The book should prove very valuable and 
stimulating for instructional and reference 
purposes in biochemistry, plant and animal 
physiology, endocrinology, pharmacology and 
medicine. Organic chemists also should find 
the book interesting. Collection of all this 
material here should help spur further work 
in this growing and useful field of chemistry. 

S. B. MAINTHIA 
The Franklin Institute Laboratories 


[HEORY OF MECHANICAL VIBRATION, by Kin 
N. Tong. 348 pages, diagrams, 6 X 9 in. 
New York, John Wiley & Sons, Inc., 1960. 
Price, $9.75. 

This book is a welcome addition to the 
extensive literature on vibration theory. It 
is essentially a condensed theory on vibration 
extending single-degree-of- 
freedom systems to multi-degree-of-freedom 
systems. Special emphasis is given to matrix 
algebra extended to generalized coordinates 
An 
interesting discussion on theory and applica- 
tion of eigenfunctions is included. 


analysis, from 


of elastic bodies and continuous media. 


The examples are well chosen, covering im- 
portant applications to transducers, lateral 
and torsional vibration of shafting, etc. In 
the analysis of a rolling load on a bridge 
beam, the author analyzes correctly the effect 


See 
: 


Oct., 1960.] 


of the inertia of the moving load in addition 
to its weight, using the interaction of the 
moving load as acting on the eigenfunction 
coordinates of the beam considered separately, 
though the author notes the equations are 
not readily solvable. 

One valuable feature of this book is the 
appendix on the matrix algebra in linear 
transformation of vectors, which use in turn 
is a characteristic feature of this book. 

The author is to be complimented on a 
remarkably complete and condensed theory 
of vibration analysis, covering a range of in- 
teresting topics. 

R. EKSERGIAN 
The Franklin Institute Laboratories 


Puysics AND MEDICINE OF THE ATMOSPHERE 
AND SPACE, edited by Otis O. Benson, Jr. 
and Hubertus Strughold. 645 pages, illus- 
trations, 6 9 in. New York, John Wiley 
& Sons, Inc., 1960. Price, $12.50. 


In November of 1951 a symposium on the 
“Physics and Medicine of the Upper Atmos- 
phere’’ was held at San Antonio, Texas. It 
was sponsored by the Air University School 


of Aviation Medicine. About a year later, 
the proceedings of this symposium was pub- 
lished with one of the editors, Brig. General 
Otis O. Benson, Jr., Commandant of the 
USAF School of Aviation Medicine. 

Seven years later in November 1958 the 
second symposium was held in San Antonio 
and the proceedings of this meeting has been 
published as Physics and Medicine of the 
Atmosphere and Space. Again, one of the 
editors is General Otis O. Benson, still head- 
ing up the USAF School of Aviation Medicine. 

With this type of editorship, we should 
expect this new book to parallel the impact 
of the first symposium proceedings. In truth, 
that is exactly what has taken place. The 
Arrangements Committee has once more in- 
vited the ranking scholars in the many dis- 
ciplines which border space flight. In the 1951 
symposium only one man was invited from 
outside the United States. The 1958 sym- 
posium saw five authorities from the free 
world present contributions. 

The symposium proceedings covered all 
facets of space flight. Some this reviewer 
considers as extraordinarily comprehensive 
and complete. As an example, Dr. S, Fred 
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Singer covers “Interplanetary Dust and Radi- 
ation Environment."’ Because of the re- 
viewer's research activities he was particularly 
interested in this section which is magnifi- 
cently treated. Dr. Leslie R. Shepherd con- 
tributes a ‘Survey of Possible Propulsion 
Systems” which again is excellent exposition. 
With characteristic thoroughness he covers 
these systems from the chemical to the photon 
rocket. The indefatigable Krafft A. Ehricke’s 
contribution on “‘Manned Orbital and Lunar 
represents a rather complete 


Space Vehicles’ 
wrap-up of man's steps into space to ulti- 
mately land on the moon. This contributor 
again stresses the need for the use of nuclear 
power plants in these space missions but it 
should be remembered that this contributor 
has for a long time advocated this approach. 

With the astronomical background of this 
reviewer it is natural for him to give particular 
attention to the papers by Dr. Roberts on 
the Sun; Dr. Kuiper on the lunar and plan- 
etary environments and Dr. de Vancouleurs 
on Mars. 

The pace of astronautics is moving forward 
so rapidly that it leaves even the experts in 
its wake. As examples of this progress, Dr 
Van Allen, in his article on radiation hazards, 
speaks of space flight at relativistic velocities 
and shows that the radiation exposure level 
reaches 20,000,000 roentgens per hour. He 
also indicates the erosive effect of particles 
on the space ship at this speed. This re- 
viewer has considered this problem in a recent 
book and concludes that the only solution to 
the problem of relativistic velocities is the 
use of an asteroid with walls several hundred 
feet thick. 

Also we find in this book in two places a 
perpetuation of a fallacy. Drs. Gerathewohl 
and Ward in their “Psychophysiologic and 
Medical Studies of Weightlessness’’ and Dr 
Stapp in his “Human Tolerance to Acceler- 
ations of Space Flight" indicate that Laika 
was weightless in Sputnik II 

This reviewer has studied this problem in 
light of the rotation of the capsule at about 
once every four seconds and assuming the 
vestibular organ of the dog was 5 feet from 
the center of rotation. He finds that the 
centrifugal force to which Laika was subject 
is the equivalent of a 4 g field! Thus, unless 
the head of Laika was precisely at the center 
of rotation she was not weightless and the 
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physiologic responses were not compatible 
with the absence of a g field. 

On the whole, this companion volume to 
“Physics and Medicine of the Upper Atmos- 
phere” is a most welcome addition for the 
serious student and worker in the space flight 
sciences. The Arrangements Committee and 
Editors are to be congratulated on the uni- 
form excellence of this volume 

I. M. Levitt 
The Fels Planetarium 
of The Franklin Insvitute 


CrassicaL Dynamics, by R. H. Atkin. 273 
pages, diagrams, 54 84 in. New York, 
John Wiley & Sons, Inc., 1959. Price, $5.25. 


Classical dynamics and electrodynamics 


constitute the essential basis of classical 
Yet the reiative time alloted them 
in academic curricula has been progressively 


reduced, in view of the rapid growth of many 


physics. 


modern branches of physics (atomic, nuclear 
and solid state physics, for example). 
the need for well written and concise presenta- 
tions. To the excellent treatises of Goldstein, 
Synge and Griffith, etc., Atkin has added a 
new one which, in a sense, is more old fash- 


Hence 


ioned and more elementary 

The vector notation is used throughout, 
and the mathematical preliminaries are intro- 
duced in the first two chapters. The central 


BOOK 


THEORY OF THERMAL STRESSES, by Bruno A 
Boley and Jerome H. Weiner. 586 pages, 
diagrams, 6 9in. New York, John Wiley 
& Sons, Inc., 1960. Price, $15.50 


This volume gives a comprehensive account 
of all the presently available theoretical tech- 
niques for determination of stresses produced 
in a solid object when it is subjected to spec- 
ified heating conditions. Thermal stress is 
explained through the basic subjects of ther 
modynamics, heat-transfer theory, elasticity 
and inelasticity. Realizing that many persons 
involved in thermal-stress analysis have little 
familiarity with heat-transfer theory, a self- 
contained review of the subject—together 
with a discussion of methods of problem form- 
ulation and methods of their solution—has 


In addition to the presenta- 


been included. 


NOTES 


force problem is presented with clarity and 
completeness, and so are the topics of rigid 
body dynamics, and normal coordinates. The 
examples solved are well chosen and illustrate 
directly the principles. The numerous prob- 
lems proposed at the end of the chapters, 
taken mostly from recent English examina- 
are rather traditional and narrow in 
They will nevertheless contribute to 
knowledge to the 


tions, 
scope. 
give a 

students. 


good working 

It is however regrettable that some topics 
of considerable importance in modern physics 
have mention the 
scattering 
Theory 
others. 


omitted. Let us 
problem, the Hamilton-Jacobi 


been 


and relativistic dynamics among 
The rich illustrations provided by 
electrodynamics have been neglected. Re- 
grettable also and rather surprising is the lack 
of an index. 

Despite these omissions, which can be rem- 
edied by an instructor, the topics considered 
are treated with an unusual skill. The work- 
ing knowledge has been emphasized at the 
expense of the philosophical foundations. As 
it stands this little book is eminently teach- 
able and could constitute an excellent text 
at the junior-senior level. 

HENRI AMAR 
The Franklin Institute Laboratories 
and Temple University 


tion of the theoretical foundations, methods 
of practical thermal-stress analysis, utilizing 
strength of materials approximations, are 
treated fully in both the elastic and inelastic 


domains. 


SELF-SATURATING MAGNETIC AMPLIFIERS, by 
Gordon E. Lynn, Thaddeus J. Pula, John 
F. Ringelman and Frederick G. Timmel. 
217 pages, diagrams, 6 9 in. New York, 
McGraw-Hill Book Co., Inc., 1960. 
$8.00. 


Price, 


This thorough treatise on self-saturating 
magnetic amplifiers, on facts of their operation 
and on precise accurate design techniques, 
shows the reader how to design complex mag- 
netic amplifier circuits with a minimum of 
Clearly ex- 


cut-and-try experimentation 
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plained in line with readily observable lab- 
oratory phenomena is the true operation of 
a physically realizable amplifier, and ideali- 
zations of nonlinear functions are made only 
after the real conditions are discussed. Among 
the subjects covered are the basic assumptions 
involving circuital laws, modern techniques 
for dealing with the transient behavior of 
self-saturating magnetic amplifiers, trans- 
former theory, feedback control theory, and 
magnetic theory. In addition, there are facts 
on the design of magnetic amplifiers, testing 
methods applicable to magnetic core material, 
and a number of commonly encountered de- 


sign problems. 


TRANSACTIONS OF THE VACUUM METALLURGY 
CONFERENCE 1959, Edited by Rointan F. 
Bunshah. 212 pages, illustrations, 8} 11 
in. New York, Associated College Presses, 
1960. Price, $7.50 (paper). 


This volume is a collection of the papers 
presented at the third annual meeting on 
vacuum metallurgy held at the University 
Heights campus of New York University, 
June 1-3, 1959, under the sponsorship of the 
Department of Metallurgical Engineering in 
co-operation with the Office of Special Services 
to Business and Industry. This conference, 
at which papers on new research and advanced 
engineering information in vacuum metallurgy 
were presented, followed in logical succession 
the two previous meetings which were courses 
consisting of a planned series of lectures 
designed to familiarize the engineer with the 
vast and diverse areas embraced by the term 
“vacuum metallurgy.” 


114 pages, 
New York, Prentice- 
(cloth); $1.50 


THe by Carl P. Swanson. 
illustrations, 6 * 9in. 
Hall, 1960. Price: $2.95 
(paper). 

First in the Foundations of Modern Biology 
Series, this text opens up the whole fascinating 
world of cellular phenomena and gives the 
reader an idea of the place of the cell as the 
basic unit of biological structure and function. 
rhe cell concept is discussed as it was orig- 
inally proposed and as it is now interpreted. 
It is considered in terms of its origin through 
cell division, correlates structure and function, 
and traces the cell through the processes of 


differentiation, reproduction and 


growth, 
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death. The Cell, along with its ten companion 
volumes in the series soon to be published, 
constitutes a new and vital contribution to 
general biology literature. This clearly written 
and well illustrated book gives an outstanding 
treatment of its own subject area, and at the 
same time, is an integral part of the whole 


series. 


ANIMAL GROWTH AND DEVELOPMENT, by 
Maurice Sussman. 114 pages, illustrations, 
6 9 in. New York, Prentice-Hall, In 
1960. Price: $2.95 (cloth); $1.50 (paper) 


This second offering in the Foundations of 
Modern Biology Series reflects the latest re- 
search findings and discoveries based on pres- 
ent-day concepts of cell physiology and ge- 
Contrasted with the usual classical 
not re- 


netics. 
treatments, the subject 
stricted merely to embryoni« 
The book begins with a consideration of single 
celled organisms and of phenomena associated 
with the development of a cell from its ap- 
pearance as a daughter cell to its disappear- 
The life cycles of typical 
microorganisms are also described as well as 


matter is 
development 


ance as a mother. 


sexual phenomena, morphogenetic processes, 
and persistent rhythmic developments are 
discussed. Readable, authoritative and ac- 
curate, this valuable study refers throughout 
to the major unsolved problems of 20th cen- 
tury biology and points out the exciting work 
which remains to be accomplished by research 
scientists of the future. 


THE CORRESPONDENCE OF ISAAC NEWTON 
VoLuME II, 1676-1687, edited by H. W 
Turnbull. 552 pages, illustrations, 8 « 11 
in. New York, Cambridge University 
Press, 1960. Price, $25.00. 


Volume II of this series sponsored by the 
Royal Society covers the years 1676 to 1687 
and includes, as in Volume I, correspondence 
not only ef Newton himself, but also of others 
(such as Collins, Oldenburg, Leibniz, Hooke, 
Halley) whose letters to each other bear a 
relation to Newton's own letters and other 
writings. The relationships among the letters 
are carefully noted, so that it becomes easy 
to trace through a series of letters, both 
prior to and following any given piece of cor- 
respondence. Those letters which were written 
in German or Latin are printed in their orig- 
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inal text, followed by English translations 
The eleven-year period covered by Volume 
I] is of particular interest, since it is climaxed 
by the publication of the Principia in 1687, 
and several of the items (notably by Halley 
and Newton) deal with events leading up to 
that publication. Volume II meets the ex- 
cellent standards set in the first volume, re- 
viewed in this JouRNAL for February, 1960 


OXIDATION—-REDUCTION POTENTIALS OF Or- 
GANIC SysTEMS, by W. Mansfield Clark. 
584 pages, diagrams, 6 X 9 in. Baltimore, 
The Williams & Wilkins Co., 1960 
$13.50. 


Price, 


Herein are presented more than 100 tables 
summarizing existing data on groups of or- 


ganic oxidation-reduction systems Along 
with critical notes and a bibliography, these 
tables form a valuable reference work. The 


remainder of the text (over half the book) is 
divided into thirteen chapters dealing with 
techniques, methods, apparatus, theories and 
definitions of many phases of organic chem- 
istry. The author points out restrictions of 
a theoretical nature, as well as limitations of 
particular experimental methods. 


My Lire AND TIMEs As A Puysicist, by I. I. 
Rabi. 55 pages, 54 x 84 in 
The Castle Press, 1960. Price, $2.75 


Pasadena, 


Iwo lectures by I. I. Rabi provide thought- 
ful readers an opportunity to take stock of 
the position of science in community life, now 
The first lecture, ““The 
makes a sound 


and for the future. 
Scientist and the Statesman,” 
plea for a more adequate understanding of 
the 


The second lecture on “Science 


science among leaders of our national 


government 


Book Notes 


4. 


and the Humanities at Mid-Century” stresses 
the importance of having nonscientists under- 
stand the scientific point of view (particularly 
students). Dr. Rabi develops his ideas by 
means of examples, both historical and auto- 
biograqhical, couched in lucid, plainly under- 
standable language. His purpose is to 
“.. , find ways to bring a greater unity to 
our culture .. he believes that ‘The 
appeal of science lies very deep in men, be- 
cause it expresses fundamentally the desire 
to explore, to understand, to know how things 
hang together, and not least, to utilize 
this knowledge for human purpose."’ His 
conclusion is that “Only by the fusion of 
science and the humanities can we hope to 
reach the wisdom appropriate to our day and 


generation.” 


DicitaL COMPUTER CONTROL ENGI- 
NEERING, by Robert Steven Ledley. 835 
pages, illustrations, 6 X 9 in. New York, 
McGraw-Hill Book Co., Inc., 1960. Price, 
$14.50 


Intended as an elementary text for seniors 


AND 


or first-year graduate students in engineering, 
this book deals with all aspects of digital 
computers from an engineering point of view. 
However, the book will prove useful in in- 
dustry as well as in the classroom. Three 
major phases are treated: (1) the over-all 
design of digital systems; (2) the logical de- 
sign of digital circuitry; and (3) the electronic 
design of digital circuits. Among the topics 
covered thoroughly are programming and 
coding, numerical analysis, Boolean algebra, 
computational methods, matrices, arithmetic 
operations, semiconductor and magnetic ele- 
ments, etc. An appendix contains several 
useful tables of formulas 
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PUBLICATIONS RECEIVED 


PuBLIC RELATIONS AND MANAGEMENT, by David Finn. 175 pages, 5 X 74 in. New York 

Reinhold Publishing Corp., 1960. Price, $4.50. 

CAREERS IN Puysics, by Alpheus W. Smith and Winston L. Hole. Revised edition, 310 pages, 
illustrations, 6 X 9 in. Columbus, Long's College Book Shop, 1960. Price, $5.95. 

New DESCRIPTION OF THORIUM SPECTRA, by Romuald Zalubas. National Bureau of Standards 
Monograph 17, 103 pages, 8 X 10} in. Washington, U. S. Government Printing Office, 
1960. Price, $0.65 (paper). 

Scientiric Basis FOR WorLp Civitization, by Leo J. Baranski. 649 pages, 6 X 9} in, 
Boston, Christopher Publishing House, 1960. Price, $7.50. 

THE SKELETON Key OF MATHEMATICS, by D. E. Littlewood. Reprinted from the first edition 
(1949), 138 pages, diagrams, 54 X 8 in. New York, Harper & Brothers, 1960. Harper 
Torchbook 525. Price $1.25 (paper). 

ue HiGHER ARITHMETIC, by H. Davenport. Reprinted from the 1952 edition, 172 pages, 
diagrams, 54 X 8 in. New York, Harper & Brothers, 1960. Harper Torchbook 526. 
Price, $1.35 (paper). 

A SHort History or Cuemistry, by J. R. Partington. Reprinted from the third edition 
(1957), 415 pages, illustrations, 54 K 8in. New York, Harper & Brothers, 1960. Harper 
Torchbook 522. Price, $1.95 (paper). 

A Direct Entry To OrGanic Cuemistry, by John Read. Reprinted from revised (1959) 
edition, 268 pages, diagrams, 54 K 8 in. New York, Harper & Brothers, 1960. Harper 
Torchbook 523. Price, $1.50 (paper). 

CHEMISTRY IN THE SERVICE OF MAN, by Alexander Findlay. Reprinted from the eighth 
(1957) edition, 318 pages, diagrams, 54 X 8 in. New York, Harper & Brothers, 1960. 
Harper Torchbook 524. Price, $1.75 (paper). 

A STRUCTURE OF SCIENCE, by Joseph H. Simons. 269 pages, 5} K 84 in. New York, Philo- 
sophical Library, 1960. Price, $4.75. 

ENERGIES OF THE UNIVERSE, by Eugene Fritz. 124 pages, diagrams, 54 K 84 in. New York, 

Philosophical Library, 1960. Price, $4.75. 


NEW BOOKS IN THE FRANKLIN INSTITUTE LIBRARY 


MATHEMATICS 


PRINCIPLES OF ANALOG CompuTaTION, by G. W. Smith. McGraw-Hill, 1959. 
CALCULUS OF FUNCTIONS OF ONE ARGUMENT, by E. J. Cogan. Prentice-Hall, 1960 
PROBABILITY, by S. Goldberg. Prentice-Hall, 1960. 

CALCULUS, WITH ANALYTIC GEOMETRY, by R. E. Johnson. Allyn and Bacon, 1960. 


PHYSICS 


EXPERIMENTS IN ELectronics, by W. H. Evans. Prentice-Hall, 1959 
ReEacToR HANpDBOOK, U.S. Atomic Energy Commission. Interscience Pub., 1960. 
ENGINEERING MECHANICS, by D. F. Gunder. Wiley, 1959. 

NuCLEAR Fue ELEMENTS, International Symposium on Nuclear Fuel Elements. Reinhold, 

1959. 

URANIUM PRopUCTION TECHNOLOGY, by C. D. Harrington. Van Nostrand, 1959. 
Non-CRrYSTALLINE Souips, by V. C. Frechette. Wiley, 1960. 
PHOTOCONDUCTIVITY OF SoLips, by R. H. Bube. Wiley, 1960. 
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Quantum Evectronics, Conference on Quantum Electronics. Columbia University Press, 
1960. 

SOME MATHEMATICAL METHODS OF Puysics, by G. H. Goertzel. McGraw-Hill, 1960. 

PRINCIPLES OF Optics, by M. Born. Pergamon Press, 1959. 

ELECTROMAGNETIC FIELDS, ENERGY AND Forces, by R. M. Fano. Wiley, 1960. 

THEORY OF MECHANICAL VIBRATION, by K. N. Tong. Wiley, 1960. 

PROBLEMS OF Low TEMPERATURE PHYSICS AND THERMODYNAMICS, International Institute of 
Refrigeration. Pergamon Press, 1959. 

NON-HOMOGENEITY IN ELASTICITY AND PLAsticity, International Union of Theoretical and 
Applied Mechanics. Pergamon Press, 1959. 

RECENT RESEARCH IN MOLECULAR Beams, by I. Estermann. Academic Press, 1959. 

THE LLoyp WiLLIAM TAYLOR MANUAL OF ADVANCED UNDERGRADUATE EXPERIMENTS IN 
Puysics, by T. B. Brown. Addison-Wesley, 1959. 

GREAT EXPERIMENTS IN Puysics, by M. H. Shamos. Holt, 1959. 


CHEMISTRY AND CHEMICAL ENGINEERING 


FATS AND O1Ls, by H. G. Kirschenbauer. Reinhold, 1960. 

CHEMICAL INSTRUMENTATION, by H. A. Strobel. Addison-Wesley, 1960. 

\N INTRODUCTION TO STATISTICAL THERMODYNAMICS, by T. L. Hitt. Addison-Wesley, 1960. 

PHOTOCHEMISTRY IN THE LIQUID AND SOLID STATES, Symposium on Photochemistry in the 
Liquid and Solid States. Wiley, 1960. 

PRINCIPLES OF Unit OPERATIONS, by A. S. Foust. Wiley, 1960. 

FAIRCHILD’s DICTIONARY OF TEXTILES, by S. S. Marks. Fairchild Publications, 1959. 

CONFERENCE OF HYPERCONJUGATION. Pergamon Press, 1959. 

TREATISE ON ANALYTICAL CHEMISTRY, by I. M. Ko_tuorr. Interscience Encyclopedia, 1959. 

CRYSTAL CHEMISTRY OF SIMPLE CoMPOUNDS OF URANIUM, THORIUM, PLUTONIUM, NEPTUNIUM, 
by S. E. Makarov. Consultants Bureau, 1959. 

METALLIZING OF PLastics, by H. Narcus. Reinhold, 1960. 

CHEMICAL ANALYSIS OF RESIN-BASED COATING MATERIALS, by C. P. Kappelmeier. Inter- 
science, 1959. 

Gtass: Its INDUSTRIAL AppLicaTions, by C. J. Phillips. Reinhold, 1960. 


ELECTRICAL ENGINEERING 


PROCEEDINGS OF THE EIA CONFERENCE ON MAINTAINABILITY OF ELECTRONIC EQUIPMENT, 
Electronic Industries Association. Engineering Publishers, 1958. 

FUNDAMENTALS OF SEMICONDUCTORS, by M. G. ScroGGie. Gernsback Library, 1960. 

DirEcT CONVERSION OF HEAT TO ELEctricity, by J. Kaye. Wiley, 1960. 

PROPERTIES, PHYSICS AND DESIGN OF SEMICONDUCTOR Devices, by J. N. Shive. Addison- 
Wesley, 1960. 

Circuit THEORY OF LINEAR Notsy Networks, by H. A. Haus. Wiley, 1959. 

SELECTED SEMICONDUCTOR CrRcUITS HANDBOOK, by S. Schwartz. Wiley, 1960. 

MopERN NETWORK ANALYsis, by F. M. Reza. McGraw-Hill, 1959. 

[TRANSFORMERS FOR THE ELECTRIC PoweR INpustTRY, by R. L. Bean. McGraw-Hill, 1959. 

ELECTRONIC CIRCUITS, SIGNALS, AND SysTEMS, by S. J. Mason. Wiley, 1960. 

POWER STATION ENGINEERING AND Economy, by B. G. Skrotski. McGraw-Hill, 1960. 

ELECTRIC ENERGY CONVERSION, by Y. Ku. Ronald Press, 1959. 

\NALYsiIs OF ELectric Circuits, by E. Brenner. McGraw-Hill, 1959. 

INSTRUMENTATION IN SCIENTIFIC RESEARCH, by K. S. Lion. McGraw-Hill, 1959. 

ENCYCLOPEDIC DICTIONARY OF ELECTRONICS AND NUCLEAR ENGINEERING, by R. I. Sarbacher. 
Prentice-Hall, 1959. 

THE THEORY OF Optimum Noise Immunity, by V. A. Kotel'nikov. McGraw-Hill, 1959. 

HANDBOOK OF ELEcTRONIC ContROL Circuits, by J. Markus. McGraw-Hill, 1959. 

MATERIALS AND TECHNIQUES FOR ELECTRON TuBEs, by W. H. Kohl. Reinhold, 1960 
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SPACE SCIENCE 


FUNDAMENTALS OF ROCKET PROPULSION, by R. E. Wiech. Reinhold, 1960 
INTRODUCTION TO ROCKET TECHNOLOGY, by V. I. Feodos’ev. Academic Press, 1959. 
GUIDE TO THE SPACE AGE, by C. W. Besserer. Prentice-Hall, 1959. 

ROCKET PROPELLANT HANDBOOK, by B. Kit. Macmillan, 1960. 


AIRCRAFT AND MISSILE DESIGN AND MAINTENANCE HANDBOOK, by C. A. Overbey. 


millan, 1960. 


ROCKET PRopuULSION, by Marcel Barrére. Elsevier, 1960. 
Space AGE Dictionary, by C. McLaughlin. Van Nostrand, 1959. 


MATERIALS SCIENCE 


INTERNAL STRESSES AND FATIGUE IN METALS, Symposium on Internal Stresses and Fatigue in 


Metals. Elsevier, 1959. 


MECHANICS AND PROPERTIES OF MATTER, second ed., by R. J. Stephenson. Wiley, 1960. 
MATERIALS SELECTION FOR Process PLANTS, by R. E. Gackenbach. Reinhold, 1960. 
NONDESTRUCTIVE TESTING HANDBOOK, by R. C. McMaster. Ronald Press, 1959. 
FERRITES, by J. Smit. Wiley, 1959. 


VacuuM METALLURGY, by R. F. Bunshah. 


Reinhold, 1958. 


YOU CAN ADVANCE SCIENCE EDUCATION 


Today, more than ever before in its 135-year history, there is vital need for 
The Franklin Institute to effectively promote education in science and technology. 
It is imperative that we meet this challenge by providing adequate educational op- 
portunities in these fields. This requires vision, objective planning, and money. 
We have more than enough of the first two requisites, but far too little of the third. 

Our programs are aimed at professional scientists and industry, as well as the lay 
public and young people seeking inspiration and guidance in choosing a career. 
The Institute’s educational programs are impressive, for they begin with students in 
the early grades of our elementary schools and continue throughout an individual's 
professional or industrial life. With more funds at our disposal, the scope and 
vigor of these activities could be greatly increased. 

The Franklin Institute is not richly endowed. It is a non-profit organization, 
depending for encouragement and support on an understanding public. Capable 
and conservative management assures wise administration of all funds. 

Your gift or bequest, large or small, will be deeply appreciated and will be used 
effectively to broaden the Institute’s educational usefulness. There is a warm 
satisfaction in giving financial support to an organization that has pioneered in, and 
is dedicated to, the advancement of science and technology through education. 

When property is transferred, title should be in the name of The Franklin 
Institute of the State of Pennsylvania for the Promotion of the Mechanic Arts. 

The Secretary of The Franklin Institute will gladly furnish you with additional 
information. Write to him at The Franklin Institute, Benjamin Franklin Parkway 
at Twentieth Street, Philadelphia 3, Pennsylvania. 
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NATIONAL BUREAU OF STANDARDS NEWS 


INTERCOMPARISON OF LABORATORY TEST METHODS 


A method for evaluating mathematically the factors involved in the vari- 
ability of interlaboratory test results has been developed by J. Mandel and 
T. W. Lashof of the National Bureau of Standards.'? This “model” for 
statistical analysis has been successfully applied to several problems involving 
both physical and chemical properties of materials. A better appraisal of test 
methods can be made by utilizing the method to distinguish between random 
and systematic sources of error. 

While the model is designed primarily for the study of a single test method 
on an interlaboratory scale, it can also be used for the interlaboratory com- 
parison of two alternative test methods. For this purpose, it is applied in 
conjunction with another statistical procedure, the sensitivity criterion,* 
involving the functional relation between the results obtained by use of the 
two alternative test methods.‘ 

The true reproducibility and therefore the usefulness of a method of meas- 
urement are unknown until the method has been evaluated on an interlabora- 
tory scale. However, the interpretation of data from interlaboratory round 
robins is not always straightforward. There has been need for a clearly 
formulated model that would relate mathematically the various factors enter- 
ing into the measurements. In recent years, considerable attention has been 
given to studies of test methods involving relatively large numbers of labora- 
tories. These studies resulted from the observation that even when well- 
established methods of measurement are used, the agreement among labora- 
tories is, in general, much less satisfactory than the agreement which can be 
achieved within a single laboratory. 

The sources of variability usually considered in such studies include : 


1. Replication error, or differences among results obtained at one labora- 
tory by one operator using the same instrument on specimens of the same 


material ; 

2. laboratory main effect, or systematic constant differences among re- 
sults obtained at different laboratories; and 

3. laboratory-material interaction, or variations among laboratories in 
the differences among results obtained for various materials. 


The present method separates the third source of variability into two com- 
ponents; one is random, and the other reflects systematic laboratory differences 

'The Measuring Process,"” by John Mandel, Technometrics, Vol. 1, p. 251 (1959). 

2**The Interlaboratory Evaluation of Testing Methods,"’ by John Mandel and T. W. 
Lashof, ASTM Bul., No. 239, p. 53 (1959) 

* The sensitivity of a method A used in the study of a property Q is defined as (AA/AQ)/oa, 
where A is the measurement obtained and a4 is the standard deviation of method A. 

*“Evaluation of Test Methods by the Sensitivity Criterion,” NBS Tech. News Bul., 
Vol. 40, p. 139 (1956 
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that cannot be represented by constant biases. To account for the random 
scatter one must consider, in addition to replication error, the effects, upon the 
measured value, of properties other than that property being measured. The 
interfering effect of such properties often varies among laboratories, causing 
additional scatter in the interlaboratory results. Another source of random 
scatter is improper calibration of equipment. In certain instruments having 
multiple scales, part of the random scatter has been shown to arise from lack 
of continuity among the scales or from variation in the ratios among the scales. 

In the absence of systematic errors, a graph of the results from one labora- 
tory plotted against those from another shows a scatter of data points about a 
straight line bisecting the angle formed by the axes (a line of 45° slope). If 
the results from one laboratory were consistently higher or lower by a con- 
stant amount than those from the other, the line would be higher or lower but 
would be parallel to the first line; that is, it would lie at 45° to either axis but 
would not pass through the origin. 

In practice, the lines obtained are not of 45° slope, and it is this non- 
parallelism or slope effect that reflects a varying systematic difference between 
the two laboratories. For example, an error in the determination of the 
normality of a reagent will cause a proportional type of error in all analyses 
employing this reagent. For any number of laboratories, the slope effect can 
be generalized by plotting the results from each laboratory against the average 
of the results from all. 

The separation of various error sources provided by this method permits 
evaluation of possible benefits derivable from reference samples. In the case 
of constant differences among laboratories, a single reference sample is suffi- 
cient. If the slopes of the lines plotted for various laboratories differ to a 
great extent, the use of two reference samples is indicated. Suppose that 
specimens of two reference materials are supplied to all participating labora- 
tories. The two materials, although necessarily similar in nature, must differ 
considerably in the value of the property under test. The values for this 
property obtained by measurement in any given laboratory for these two 
materials can be compared with the known correct (or assigned) values; then 
at each laboratory the measured values for any submitted unknown sample 
of such material can be “adjusted”’ and brought into line with the similarly 
adjusted values from other laboratories. 

The use of reference samples to improve interlaboratory precision is an 
accepted procedure and differs little from some calibration methods; however, 
reference samples should be used for this purpose only after all economically 
feasible improvements in the control of environmental conditions or other de- 
sirable modifications of the test procedure have been exhausted. 

Most statistical methods are highly sensitive to extremely discrepant 
values. The present method contains several safeguards against misleading 
results due to the inclusion of outlying measurements. Standard deviation 
of replication error is examined for each material and laboratory. Single 
discrepant values are then easily detected from inflation of the corresponding 
standard deviations. Trends are handled statistically through appropriate 
transformations of scale. If all the replicates obtained by one laboratory for a 
particular material are affected by the same error, the contribution of this 
laboratory to the random part of the laboratory-material interaction will be 
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unduly large and consequently obvious. A laboratory discrepant in all values 
will be apparent from the fact that its line will differ from the others appreciably 
in position or slope. 

Sometimes these discrepancies can be clarified and the nonconforming 
values replaced by the results of measurements taken under proper conditions. 
When no reason for the conflict can be found, the values must be considered 
as possibilities in the application of the test method. 


A CARBON-14 STANDARD FOR LIQUID-SCINTILLATION COUNTERS 


The National Bureau of Standards has prepared a new carbon-14 standard 
sample to meet the need created by the increased use of liquid-scintillation 
counters, particularly in biochemical and medical studies. The solvents most 
often used in these counters are toluene and xylene; the standard used in their 
calibration must therefore be soluble in and compatible with these compounds. 
To fulfill this requirement, W. F. Marlow and R. W. Medlock of the Bureau's 
radioactivity laboratory used benzoic acid-7-C™ dissolved in toluene for the 
new standard.°® 

In the past few years, liquid-scintillation counting systems have become 
increasingly accurate and reliable, while their cost has greatly decreased. 
The sample of material to be counted by this method is more easily prepared 
than for most other methods, and since it is dissolved directly in the solution, 
the counting efficiency is usually much greater, especially for low-energy radia- 
tions of carbon-14. 

The sample of benzoic acid-7-C in toluene was standardized by quantita- 
tive oxidation of the benzoic acid and toluene. The carbon dioxide produced 
was collected and the level of radioactivity of this carbon dioxide was com- 
pared, in an ionization chamber, with the radioactivity of carbon dioxide 
prepared quantitatively from the national sodium carbonate-C"™ standard. 

Previous work® has shown that C™O, and CO, may be evolved at different 
rates during oxidation of organic compounds either by wet combustion or by 
burning in an oxygen stream. The possibility of error due to this isotope 
effect was avoided by burning the sample in a modified Paar oxygen bomb to 
obtain complete oxidation. 

The new carbon-14 standard sample may be obtained from the Radio- 
activity Section, National Bureau of Standards, Washington 25, D. C. for 
$27.00.? 


Toluene, for Liquid-Scintillation Counters,’’ by W. F. Marlow and R. W. Medlock, J. Research 
NBS, Vol 64A, p. 143 (1960). 

* “Errors of Combustion of Compounds for C'* Analysis,” by W. D. Armstrong, L. Singer, 
S. H. Zbarsky and B. Dunshee, Science, Vol. 112, p. 531 (1950 

7A complete listing of standard samples is contained in Standard Materials, NBS Circular 
552 (third edition), which may be obtained by writing to the Superintendent of Documents, 
U.S. Government Printing Office, Washington 25, D. C. (35 cents). 
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Satellites Inflated with Plastic 
Foam?—Tomorrow’s space _ stations 
and special-purpose satellites may be 
built with balloon-like tubes that can 
be folded into compact packages for 
launching and inflated automatically 
in space with a plastic foam hardened 
by, the sun’s heat, according to an en- 
gineer of the Radio Corporation of 
America. 

This technique is now under study 
at RCA’s Astro-Electronics Division, 
Princeton, N. J., as a possible solution 
to the problem of setting up in space 
very large structures, such as exten- 
sive solar energy collectors and radar 
antennas 60 ft. or more in diameter. 
It was described by Carl C. Osgood, 
Manager, Design Engineering, of the 
RCA division, in a paper prepared for 
presentation at the annual West Coast 
meeting of the American Astronautical 
Society. 

Mr. Osgood said the use of inflat- 
able tubes promises to overcome the 
size and weight limitations encoun- 
tered in launching folded structures 
designed to be set up in space. As 
an example, he pointed out that an 
inflatable tube structure weighing only 
32 Ib. could provide a disk 44 ft. in 
diameter capable of carrying a solar 
cell array weighing 420 Ib. 

“For smaller short-life applications, 
a gas bottle large enough to inflate 
and maintain the leakage rate of the 
array is considered satisfactory,” said 
Mr. Osgood. “But for the larger, 
longer-life devices, the maintenance of 
shape and rigidity by continuous gas 
pressure is inadmissible because of the 
weight of tankage. An arrangement 


of rigid ribs and hingelock combina- 
tions appears feasible, but is elimi- 
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nated because of the high weight of the 
many pieces of hardware and frequent 
difficulty in meeting the folding 
requirements. 

“Additional search for methods of 
inflation leads to the use of ‘foamed- 
in-place’ structural members. This 
process makes use of the increase in 
volume—from three to fifty times be- 
tween liquid reactants and the solid 
foam product. Before curing, the 
foam has some fluid properties, and 
its pressure, plus that of the escaping 
by-products, can be used as the infla- 
tion pressure. The cured foam could 
provide the permanent rigidity re- 
quired to retain the shape of the array. 
This would obviate the need for con- 
tinuous gas pressure or for transport- 
ing rigid structural members.” 

Mr. Osgood said that engineers at 
RCA’s Princeton Space Center are 
now investigating methods of foam 
inflation and the development of plas- 
tic foam materials sufficiently resis- 
tant to the radiation encountered in 
space. He discussed several possible 
methods of foam inflation based on 
combination of the reacting elements 
by their simultaneous injection into 
the inflatable structure when the sat- 
ellite is released in orbit from the 
launching vehicle. 

A further method under study, ac- 
cording to Mr. Osgood, is the use of 
“‘frozen’’ materials that would be com- 
bined in the collapsed satellite struc- 
ture before launch and kept below the 
reaction point by refrigeration tech- 
niques until the orbit is achieved. In 
orbit, solar heat would start the re- 
action, generating the foam to expand 
the structure, and curing the foam to 
the desired hardness. 
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Polaris Missile Trainer. [he first 
production-model Polaris Fire Control 
trainer has been delivered to the 
United States Navy by General Elec- 
tric Company. Designed to train 
crews aboard fleet ballistic missile sub- 
marines, the electronic device permits 
simulated submerged launchings of the 
intermediate-range Polaris. 

It is scheduled for installation 
aboard the nuclear-powered submar- 
ine Theodore Roosevelt, now under con- 
struction at Mare Island Navy Yard, 
Calif. Such trainers will be standard 
equipment aboard George Washington 
(SSB(N)598) and Ethan Allen (SSB- 
(1N)608) class fleet ballistic missile sub- 
marines, according to the Navy. 

The computer-like trainer can be 
connected to the fire-control system 
of the submarine, allowing its attack 
party to practice complete launching 
operation on the same equipment that 
would be used to fire an actual missile. 

As the crew operates the fire-control 
system, the device is so versatile that 
it creates the conditions of an actual 
missile launching. It can even simu- 
late a wide variety of malfunction for 
the attack party to cope with. 

General Electric conducted a full- 
scale team effort among several of its 
departments to design and produce 
the training device for the Naval 
Training Devices Center, Port Wash- 
ington, N. Y. 

Its General Engineering Labora- 
tory, Schenectady, N. Y., developed 
the complex electrical design at the 
direction of the G. E. Ordnance De- 
partment, Pittsfield, Mass. 

The company’s Specialty Control 
Department, Waynesboro, Va., was 
selected as fabricator of the equipment 
because of its wide range of manu- 
facturing capabilities and techniques 
and its experience in building similar 
special and defense control equipment, 

Contracts for production-model 
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Torics 
trainers have been let to the G. E., 
Ordnance Department and will be 
manufactured by the Specialty Con- 
trol Department. 


Subminiature Accelerometer.—The 
development of a new subminiature 
piezoelectric accelerometer, Model 
606, designed for use in special shock 
and vibration test applications where 
small size and light weight are the 
prime considerations, has been an- 
nounced by Columbia _ Research 
Laboratories, Woodlyne, Pa. 

The model 606 is a stainless steel 
accelerometer only } in. high, weighing 
just 1.0 gram, permitting many in- 
flight tests that could not normally 
have been conducted with standard 
size or miniature size pickups. The unit 
has a sensitivity of 2.5 pk-mv/pk-g 
and a frequency response of 5 cps to 
10 ke with its resonant frequency at 
135 ke. Maximum acceleration is 
given at 40,000 g, amplitude linearity 


+2°%, and temperature range —65° 
F. to +300° F. for standard units and 
— 65° F. to +500° F. for high temper- 
ature units, with less than +10 per 
cent variation in sensitivity through- 
out both temperature ranges. 


Satellite Observatories.—-A project 
to launch satellite observatories ca- 
pable of taking very accurate astronom- 
ical measurements from high altitudes 
in space is under way at the University 
of Wisconsin, one of five institutions 
participating in NASA's orbiting ob- 
servatories program. ‘The others are 
Harvard, Princeton, the Smithsonian 
Institution, and NASA’s own Goddard 
Space Laboratories. For the immed- 
iate future, in the cooperative effort, 
work will be concentrated on launch- 
ing three satellites to carry telescopes 
of various sizes and to be devoted to 
problems concerned with ultra-violet 
light. 
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At altitudes in excess of 100 miles, 
the ‘‘flying”’ observatories will escape 
the problems presented by the Earth's 
atmosphere. They should be able to 
see the ultra-violet light that comes 
from the stars and to transmit their 
measurements to earth by electronic 
methods. The telescopes should be 
able to point at a star with more ac- 
curacy and precision than earthbound 
ones. 

The first of the orbiting observa- 
tories will carry a number of small 
telescopes for UW experiments on the 
brightness of stars in the ultra-violet. 
Wisconsin will share this satellite with 
the Smithsonian Institution which will 
be doing photography of the skies by 
television techniques in the ultra- 
violet. The second satellite will carry 
an ultra-violet spectrograph for the 
Goddard Space Labs; the third, a high 
dispersion ultra-violet spectrograph 
for Princeton. 

It is planned to lift the observatories 
into orbit with Atlas rockets topped 
with Agena B boosters. A number 
of these Agenas have already been 
successfully launched in other space 
experimentation. Wisconsin hopes to 
get the first of the orbiting observa- 
tories into space by 1963. UW scien- 
tists hope, also, that once orbited, a 
satellite will continue to operate sat- 
isfactorily for at least a year. 

UW astronomers will be responsible 
for the design of the telescope instru- 
mentation and for everything neces- 
sary to carry out the actual experi- 
ments. They will construct experi- 
mentally some parts in order to test 
out certain ideas. 

Though the first satellite will carry 
relatively small telescopes for Wiscon- 
sin investigators, UW astronomers ex- 
pect to launch later in the program 
a single telescope as much as 36-in. 
in diameter and something in the 
order of a ton in weight. 
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Ever since the beginning of astron- 
omy, scientists have had to observe 
from the bottom of a ‘‘sea”’ of air and 
this atmosphere, more than any other 
thing, has limited what can be learned 
about the nature of the universe. It 
bends the light, making the images of 
the stars being measured fuzzy and 
seeming to move. Measurements in 
the finest detail are not possible. The 
blanket of air also scatters the light, 
thus producing a sky brighter than 
it would be otherwise. Because of 
this, the faintest objects cannot be 
seen. In addition, the atmosphere 
presents the problem of being opaque 
to most electro-magnetic radiation. 

It is possible to escape the bending 
and scattering of light by making ob- 
servations from balloons, but only by 
reaching altitudes in excess of 100 
miles (by means of the flying obser- 
vatories) can astronomers see the ultra- 
violet that comes from the stars. So, 
the first programs in space astronomy 
will be devoted to measuring ultra- 
violet and radio light. 

Possibly within the next year, NASA 
will precede the actual launching of 
the orbiting observatories with astro- 
nomical measurements taken in rela- 
tively short periods of observation 
from rockets or from an X-15, manned 
rocket plane. By the fall of 1961, the 
X-15 is expected to reach the heights 
necessary. Whatever the vehicles, 
rockets or planes, most of them will 
carry photo-electric equipment and 
the records they gather will be radioed 
to Earth. 

All of the institutions concerned 


with the orbiting observatories will 
ultimately share the use of the ‘‘flying”’ 
telescopes, and astronomers not as- 
sociated with the project will also have 
opportunities to make observations 
with the satellites. 
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Test-Pilot Monitor.—A new manned- 
vehicle monitor that can tell at a 
glance the immediate condition of a 
pilot and craft he is flying has been 
designed for the Air Research and 
Development Command (ARDC) by 
International Telephone and Tele- 
graph Laboratories, San Fernando 
Valley, Calif. 

The new device, to be operated by 
ground personnel, for the first time 
will give direct read-out or “‘real time” 
safety reports on the pilot’s physio- 
logical, and environmental condition. 
In addition, “how goes informa- 
tion on the vehicle itself is transmitted 
to the ground observers. 

Designed and in development under 
contract for ARDC’s Air Force Flight 
Test Center, Edwards AFB, Calif., 
the space-data display unit shows the 
changing environmental and physio- 
logical conditions of a man before he 
may be aware of them. Safety con- 
trol from the monitoring point will be 
instantaneous by communication to 
the vehicle. 

According to Paul Sehnert, AFFTC, 
Project Engineer, and Carlos Miller, 
Executive Engineer of ITTL, the 
“real time’ console will keep abreast 
of man’s conditions within fractions 
of a second and can indicate immedi- 
ately when he is undergoing or enter- 
ing an unsafe condition. No compu- 
tation of data, graph, or gauge anal- 
ysis is necessary, they said. 

“We need to know immediately the 
status of a man’s condition,”’ Sehnert 
said, ‘‘and with this new device com- 
plex information is shown and read 
without need of lengthy and perhaps 
very costly analysis.” 

The development work currently 
under way follows a feasibility study. 
The system is specifically being built 
to monitor 


experimental high-per- 


formance aircraft or rocket 


Hights such as the X-15. 


plane 
It may be 
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used on other space-type vehicles or 
man-in-space programs. 

The electronic design of the new 
system follows the most advanced 
state-of-the-art. Solid state compo- 
nents are used wherever practicable, 
and miniaturization and reliability are 
exercised throughout the console. 

The first unit will be ready for in- 
stallation and operation early next 
year at Edwards AFB. 

The console is separated into two 
units, vehicular and_ physiological, 
each about half the size of an office 
desk. A flight test engineer will mon- 
itor the vehicular unit while a flight 
surgeon will monitor the physiological 
unit, but they will only indirectly, by 
voice communication, exercise control 
of the vehicle through advice to the 
pilot. 

The information display on the ve- 
hicular panel includes various stability 
and aerodynamic parameters, struc- 
tural strains, flutter, power plant be- 
havior, and others such as skin tem- 
perature, fuel, fire warning, electrical 
equipment, etc. 

The information shown on the phys- 
iological panel includes the man’s 
heart rate, blood pressure, tempera- 
ture, breathing rate and volume, the 
parameters concerning his environ- 
ment, oxygen supply, acceleration, 
radiation count, etc. 

Four basic kinds of display media 
are used to report the so-called “how 
goes it’’ information: Cathode ray 
tubes generate, for two inter-related 
parameters, a spot of light on an X-Y 
plot, and for a single parameter, a 
marker line on a thermometer-type 
display ; indicator lights show ‘“‘safe”’ 
or “unsafe” conditions of certain pa- 
rameters ; counter wheels show quanti- 
tative value of selected parameters; 
and time-history displays, for certain 
physiological parameters, show recent 
history to allow accurate interpreta- 
tion of the real-time value. 
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The cathode ray displays will pre- 
sent a green color indication of the 
“OK” condition which changes to a 
harsh orange when either the actual 
value or the rate-of-change becomes 
unsafe. For the X-Y plot displays of 
inter-related parameters, the green 
spots of light are located inside “‘safety 
envelope” areas which have been pre- 
viously determined for the particular 
parameters for both the man and the 
vehicle. The spots of light change 
their color to orange when the pa- 
rameters move outside the ‘‘safety en- 
velope”’ areas. 


Rotating Microscope Stage.—Said 
to have almost unlimited possibilities 
in laboratory and _ production-line 
work, a brand-new Rotating Micro- 
scope Stage now permits microscopic 
top-and-side examination of tiny sub- 
jects without moving them or holding 
by hand. Developed by Edmund Sci- 
entific Company for use with low power 
microscopes, the new Rotating Stage 
is now finding use for the microscopic 
inspection and study of tiny compo- 
nents and assemblies, material sam- 
ples, contact lenses, resistors, transis- 
tors, miniature gears, gems, minerals, 
biological and other specimens. Sub- 
ject is placed on the revolving stage 
for study under the scope. A touch 
of the finger revolves the center plat- 
form and rotates the subject for ex- 
amination. An angled first surface 
mirror permits top-and-side viewing 
without the necessity of ‘‘upsetting”’ 
the subject or awkward holding by 
hand. Flat edge of rotating wheel in- 
dicates complete revolution. Unit is 
sturdy metal with plastic rotating 
parts—size 2 in. X 2 in. 2 in. 


World’s Fastest Switching Device. 
—The world’s fastest switching de- 
vice—a_ radically improved tunnel 
diode that could conceivably count ten 
billion one-dollar bills in a second or 
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make computer ‘‘decisions"’ with equal 
speed—has been introduced by the 
Radio Corporation of America. 

Already, the Department of Defense 
is studying these new tunnel diodes 
with a view to developing a computer 
that will operate at a speed of a billion 
cycles per second—100 times faster 
than any computer in existence, ac- 
cording to Kenneth M. McLaugh- 
lin, Manager, Computer Products, 
RCA Semiconductor and Materials 
Division. 

“The device is one of four new types 
already in production and available 
commercially,’ he said. ‘“‘Because of 
their tremendous speed, they are ex- 
pected to have wide use in a variety 
of electronic systems including data 
processing, missiles, satellites and var- 
ious communications equipment.” 

Tunnel diodes—tiny solid-state de- 
vices no bigger than the head of a 
match—control the flow of electrons 
that make up an electrical current 
much faster than an electron tube or 
a transistor. They are called tunnel 
diodes because of the manner in which 
electrons seem to “tunnel” through 
them with the speed of light. 

“These new diodes operate at in- 
credible speeds,”” Mr. McLaughlin 
said. “One of them is fifty times 
faster than the fastest transistor. As 
an illustration of their speed, Mr. Mc- 
Laughlin said the fastest of the new 
diodes can switch from one voltage 
to another in the time it takes light, 
traveling at 186,000 miles per second, 
to traverse a distance of about one 
inch, 

“Theoretically, then, it is possible 
for this device to switch—or, in pop- 
ular terminology, ‘to make decisions’ 
—ten billion times a second.”’ 

The low prices listed for the new 
diodes are a result of both technolog- 
ical advances that have reduced pro- 
duction costs and a growing market 
resulting from an awakening to the 


4 


CURRENT 


tremendous versatility, reliability and 
promise of these rugged and very tiny 
components. 


Repair Stations for Spacecraft. 
Emergency drydocks and ‘space 
hogs” have been proposed for fixing 
spacecraft “‘flat tires.”” The problems 
of fixing an ailing spacecraft and pos- 
sible solutions were discussed in a 
technical paper by E. J. Merrick, a 
project engineer in the General Elec- 
tric Company’s Missile and Space Ve- 
hicle Department at Philadelphia. 
The paper was presented by Dr. P. 
R. Knaff, a Department engineering 
psychologist, at the annual West 
Coast meeting of the American Astro- 
nautical Society. 

“The drydocks could be as simple 
as a plastic ‘sausage skin’ drawn over 
and around the craft and then inflated, 
or they could be as complex as a geo- 
desic sphere,’’ Merrick said. 

The drydocks would be inflated 
with gases fed from the mother ship 
or from an independent source, pro- 
viding an atmosphere for sustaining 
life. When there is work to be per- 
formed outside of the spacecraft, the 
mother ship would be brought into 
thedrydock. Light comfortably suited 
men could then work inside the sealed 
“sausage.” 

After completion of the tasks, the 
drydock would be pumped out by 
mechanical collapsing of the structure 
for return aboard the mother ship. 
The atmospheric gases of the drydocks 
would also be returned to the mother 
ship for purification. 

“Safety, reliability, comfort, and 
efficiency can be increased through 
the use of such space shelters,"’ Mer- 
rick said. 

He explained that the “space hog” 
is a General Electric concept which is 
designed to provide a comfortable, 
simple, and highly useful means for 
excursions outside of the spacecraft. 
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Its name is derived from the cassions 
used by sandhogs. 

The device consists of a pressure- 
tight, non-rigid tube or cylinder car- 
rying stiffening rings at selected in- 
tervals. When not in use, the device 
is completely stored in an airlock on 
the spacecraft. 

When the device is to be used, it 
is extended and pressurized. Then 
the spaceman, wearing a light-duty, 
comfortable suit, crawls into the space- 
hog, which can be rotated and/or 
flexed so that the working end can 
sweep over any portion of the mother 
ship to permit close-on, safe task per- 
formance. The outboard end of the 
spacehog would carry such devices as 
a viewing port, lights, tools, mechani- 
cal arms, or an air-lock with traction 
attachments. 


Molecular Biophysics Lectures.- 
To inquire into the possible physio- 
logical role of fast transfer reactions in 
eliciting specific interactions in ordered 
macromolecular structures, a_ series 
of seminar lectures was held at the 
Massachusetts Institute of Technol- 
ogy during the spring term of 1960. 
Included were lectures on electron and 
energy transfer mechanisms (Kasha, 
Taube, Murrell, Weber, Stryer) and 
on current concepts of proton trans- 
port mechanisms (Eigen and Grun- 
wald) of immediate biophysical and 
biochemical significance. Properties 
of water and ice crucial for proton 
transfer were reviewed (Frank, Brady, 
Klotz, Six lec- 
tures of a more general nature dealt 


Fernandez- Moran). 


with aspects of macromolecular inter- 


action properties (Gergely, Berendsen, 
Oncley, Davison, Wiener, Elsasser). 

Abstracts of these lectures and per- 
tinent references are available gratis, 
from Professor Francis O. Schmitt, 
Department of Biology, M.I.T., Cam- 
bridge 39, Mass. 
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INNER-SPACE KNOWLEDGE FOR OUTER-SPACE TECHNOLOGY ADVANCED BY TI'S 


scientific key to Davy Jones’ locker! 


At a time when the focus is on outer space — a time when 
we know more about the moon’s surface than we know 
about 70% of the earth’s surface — physicists such as 
this sea-going scientist from Texas Instruments are fur- 
— space-age technology with studies of inner space 
— the oceans. 


He represents one of more than 60 global TI geophysical 
investigation parties who are applying 30 years’ expe- 
rience in earth sciences to measurement of earth’s land 
and sea physical characteristics. Fundamental to modern 
defense, earth sciences have extended their traditional role 
in petroleum exploration to gathering information vital 
to programs in underwater warfare, missile accuracy, 
seismic communications, detection/surveillance of nuclear 
testa, and design of underground defense structures. 


TI's talent for these “down-to-earth” studies dates back 
to 1930, when the company was formed as Geophysical 
Service Inc. — the first independent company to perform 
seismic reflection surveys for oil exploration. And from 
GSI’s early need for sensitive electronic seismic equipment 


TEXAS 


came an engineering skill that was applied to Anti-Sub- 
marine Warfare systems in World War II, and resulted in 
TI's Apparatus division becoming a leading ASW equip- 
ment manufacturer today. This rare blending of earth 
sciences knowledge and experience with skills , # systems 
and components gives TI a unique capability in space-age 
technology. 


The founding company (GSI) is now the geophysical 
exploration arm of TI's Geosciences & Instrumentation 
division, and these capabilities are extended to govern- 
mental agencies by the Geosciences department. Designing 
and manufacturing advanced instruments and systems in 
wide use by these parties and others is the Instrumenta- 
tion Product group, with such products as WORDEN* 
gravity meters, EXPLORER® seismograph systems and 
seisaMAC® seismic computers. 


For more about geosciences at TI, write to Central Mer- 
chandising for Bulletin DM-101. 


*Trademark of Texas Instrumente 
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3 floor from charges exploded on land or in water as much 
miles away. Similar land-based 1! geophysical parties are now = 
3 working in more than 22 countries of the world 
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MATH 
ENGINEERING 


Now in the fourth year of operation, The Franklin Institute Com- 
puting Center has shared its technical know-how, modern electronic 
equipment and experienced personnel with hundreds of progressive 
industries and government agencies across the nation. 


We offer the services of creative people, highly skilled in their respec- 
tive fields and ably trained in the application of these skills to the ever 
expanding area of electronic computers and data processing systems. 


This staff is now available for analysis, system design, programming 
or coding of projects of unlimited scope or context. Input to our large 
scale computer and peripheral equipment is acceptable in any form. 
Results are provided on cards, plastic or metallic tape, and in com- 
pletely edited printed copy. Our extensive library of automatic coding, 
engineering, data processing and mathematical routines is available to 
all users, and machine time is provided with or without the services 
of programming personnel. 


THE FRANKLIN INSTITUTE 
Computing Center 


20th Street and Benjamin Franklin Parkway 
Philadelphia 3, Pa. LOcust 4-3600 
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